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PREFACE 


This  is  a  final  report  for  a  project  entitled 
"Heat  Transfer  from  an  Ionized  Gas  to  a  Gaseous  Coolant". 

Work  was  sponsored  by  the  Advanced  Research  Projects  Agency 
and  performed  under  contract  with  the  Office  of  Naval  Research, 
(Contract  No.  Nonr  1858(31))  during  the  period  February  1959 
to  30  September  1962.  During  this  period  two  semi-annual 

It 

report s(l),  (2)  ,  nine  quarterly  letter  reports,  and  four 
technical  reports  (3),  (4),  (5),  (6),  were  issued.  This  report 
summarizes  the  material  presented  in  these  earlier  reports 
and  generally  presents  the  status  of  the  project  as  of 
30  September  1962. 

The  project  has  continued  under  the  sponsorship  of 
the  Aeronautical  Research  Laboratories,  Wright-Patterson  Air 
Force  Base,  Contract  No.  AF33(657)9962.  Some  of  the  research 
performed  in  the  continuation  study,  included  in  this 
report,  were  pertinent  to  the  clarification  of  the  earlier 
work. 
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Numbers  in  parenthesis  indicate  references  at  the  end  of 
the  report. 
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I.  SUMMARY 

Theoretical  and  experimental  studies  were  made 
of  subsonic  mixing  between  a  partially-ionized  gas  jet  and 
a  cool,  coaxial  gas  flow.  The  gases  studied  were  argon 
and  helium.  Both  turbulent  and  laminar  mixing  regimes 
were  investigated.  In  the  turbulent  mixing  study,  analytical 
expressions  were  closely  verified  by  experiment.  In  laminar 
mixing,  theoretical  transport  properties  were  computed  and 
used  to  analyze  the  laminar  mixing  case. 

In  the  turbulent  mixing  experiments  argon  was 
heated  to  temperatures  up  to  14,000°K  by  a  commercial  arcjet 
and  exhausted  as  a  3/4-inch  diameter  jet  into  a  coaxial  flow 
of  helium  at  room  temperature.  A  water-cooled  calorimetric 
probe  1/8-inch  in  diameter  was  developed  and  used  to  survey 
the  jet  core  and  mixing  regions  for  enthalpy,  gas  composition 
and  stagnation  pressure.  Temperature,  velocity  and  composition 
profiles  of  the  jet  and  mixing  region  were  obtained  with 
good  accuracy,  as  indicated  by  detailed  surveys  providing 
mass  and  energy  balances  of  the  jet.  Good  agreement  between 
measured  data  and  analytical  studies  was  obtained,  with  the 
analysis  correctly  predicting  the  spreading  boundaries  of  the 
jet  and  the  axial  decay  of  energy,  momentum,  and  argon 
concentration . 

In  support  of  both  the  laminar  and  turbulent  studies 
a  detailed  analysis  of  the  approach  to  equilibrium  of  the 
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jet  was  performed,  and  the  radiation  from  the  jet  was 
studied.  It  was  concluded  that  electrons  were  close  to 
thermal  equilibrium  with  other  components  of  the  jet  in  the 
zone  studied.  An  experimental  determination  of  the 
radiation  loss  confirmed  theoretical  estimates  that  between 
5  and  10  per  cent  of  the  jet  energy  was  lost  by  radiation 
in  the  temperature  range  of  interest. 
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II.  INTRODUCTION 

The  general  purpose  of  this  work  was  to 
establish  methods  of  dealing  with  the  very  high 
temperature  gas-gas  interactions  required  by  certain 
advanced  space  propulsion  concepts.  Specif ically, the  effects 
on  conventional  gas  mixing  and  heat  transfer  theories  of 
(a)  high  ionized  fraction,  (b)  high  temperature  gradients,  (c) 
electron  equilibrium  and  (d)  radiation  from  the  jet  were 
investigated. 

The  work  was  performed  to  yield  an  understanding 
of  the  basic  transport  phenomena  and  not  with  orientation 
to  any  immediate  application.  Nevertheless,  it  is  expected 
that  work  accomplished  in  this  study  and  in  similar  studies 
will  have  eventual  application  in  several  advanced  propulsion 
systems. 

One  of  the  possible  applications  is  the  gaseous- 
core  nuclear  rocket  (e.g.,  see  References(7) ,  (8),  (9)),  in 
which  it  is  necessary  to  transfer  energy  from  a  heavily- 
ionized  .fissionable  fuel  plasma  to  a  propellant  gas  such  as 
hydrogen  or  helium.  Economic  considerations  require  that 
only  very  small  leakage  of  fuel  ions  be  allowed  (7).  Hence 
it  is  essential  that  all  the  nuclear  energy  developed  be 
delivered  to  the  light  propellant  gas,  the  ionized  fuel 
being  retained  by  magnetic  fields,  centrifugal  separation, 
or  some  combination  of  these  effects.  Performance  of  these 
powerplants  is  therefore  dictated  by  the  energy-exchange 
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processes  between  the  ionized  fuel  atoms  and  the  propellant 
gas. 

A  second  application  occurs  in  case  it  becomes 
necessary  to  obtain  short-term  high  thrust  levels  from  a 
low-thrust,  high-specific  impulse  propulsion  system  (e.g., 
thermal  arcjet).  This  can  most  easily  be  done  by  diluting 
the  high-energy  jet  with  additional  cold  propellant,  thus 
reducing  the  effective  exhaust  velocity  but  increasing  the 
thrust.  Here  the  attainable  performance  is  a  critical 
function  of  the  mixing  process  between  cold  and  hot  gases. 

A  third  application,  which  may  turn  out  to  be 
the  most  important,  is  the  film  cooling  of  chambers  for 
extreme-energy  propulsion  systems.  These  would  be  either 
of  the  gaseous-core  nuclear  type  discussed  earlier,  or, 
possibly,  magnetohydrodynamic  accelerators  using  high-temperature 
ionized  gases. 

In  general,  it  is  often  necessary  to  know  the 
characteristics  of  extreme-temperature  gas  flows.  Rocket 
exhaust  nozzles  at  the  low  pressures  required  for  space 
operations,  arcjet  engine  exhausts,  low-pressure  flows 
through  heat-transfer-type  nuclear  reactor  cores  slated 
for  space  propulsion  systems,  etc.  are  examples  of  non¬ 
equilibrium  flows  whose  analysis  is  not  yet  possible  because 
of  a  lack  of  high-temperature  gas  data. 
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III,  THEORY  AND  ANALYSIS 
A.  Turbulent  Mixing 

1.  Introduction 
History 

After  the  initial  work  of  0.  Reynolds  in  1883 
on  the  laws  of  instability  of  streamline  motion  and  the 
basic  ideas  of  turbulent  flow,  the  first  major  study  of 
turbulence  was  due  to  G.  I.  Taylor  (10)  in  his  1920  paper 
"Diffusion  by  Continuous  Movement".  However,  this  theory 
treats  only  the  case  of  isotropic  turbulence  and  cannot  be 
used  on  even  relatively  simple  flow  problems  with  any  great 
degree  of  success. 

In  contrast  to  Taylor's  statistical  approach, 
certain  phenomenological  theories  were  later  developed,  the 
best  known  of  these  being  the  mixing  length  theory  due  to 
L.  Prandtl  (11).  Tollmeim  (12),  in  1926,  applied  the  Prandtl 
mixing  length  theory  to  the  following  incompressible  cases: 
(a)  The  mixing  of  a  parallel  stream  with  an  adjacent  fluid 
at  rest,  (b)  The  mixing  of  a  two-dimensional  jet,  issuing 
from  a  very  narrow  opening,  with  a  medium  at  rest. 

In  1941  Reichardt  (13)  introduced  the  "constant 
exchange-coefficient"  theory  as  an  improvement  upon  the 
Prandtl  "mixing  length"  theory.  This  essentially  corresponds 
to  the  assumption  of  a  constant  exchange  coefficient  over  the 
entire  mixing  region  of  a  free  jet.  All  the  above  work 
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considered  the  case  of  incompressible  turbulent  mixing 
of  a  single  medium,  and  was  therefore  concerned  only  with 
momentum  transfer. 

Experimental  work  by  Forstall  and  Shapiro  (14) 
on  the  mixing  of  a  free  jet  with  an  external  gas  of  a 
different  type  then  showed  that  the  transport  of  specie 
concentration  occurs  more  rapidly  than  does  the  transport 
of  momentum.  W.  Warren  (15)  in  1957  further  extended  the 
study  of  free  turbulent  mixing,  using  an  integral  analysis 
of  a  compressible,  heated  jet  issuing  into  a  medium  at  rest. 

As  a  result  of  comments  by  Prandtl  (16)  on  the  Reichardt  "constan 
exchange  coefficient"  theory,  a  "modified  exchange-coefficient" 
concept  was  employed.  This  concept  requires  the  exchange 
coefficient  to  be  constant  across  a  given  axial  station 
but  to  vary  with  axial  position.  Agreement  with  the  data  of 
Corrsin  and  Uberoi  (17)  for  axial  temperature  decay,  jet 
temperature  spreading,  axial  velocity  decay,  and  jet  velocity 
spreading  was  excellent.  However,  the  temperature  range 
was  quite  small  ,  extending  only  about  500°R  from  the  room 
temperature  reference.  This  study  was  done  for  subsonic 
and  supersonic  jets  with  both  primary  and  secondary  gases 
of  the  same  species.  Velocity  and  temperature  profiles  were 
also  measured  at  various  axial  positions. 

Existing  Phenomenological  Theories 

The  study  of  turbulent  flow  problems  generally 
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requires  the  application  of  one  of  the  following  semi- 
empirical  theories: 

1.  Momentum  transfer  theory  due  to  Prandtl. 

2.  Modified  vorticity  transfer  theory  due  to 
Taylor. 

3.  One  of  the  various  exchange  coefficient 
theories  due  primarily  to  Reichardt. 

The  basic  purpose  of  each  of  these  theories  is 
to  relate  the  turbulent  shear  mechanism  to  mean  flow 
properties  so  that  the  differential  or  integral  equations 
describing  the  problem  may  be  solved.  In  a  sense  each  of 
these  theories  is  somewhat  unrealistic  (i.e.,  it  has  been 
experimentally  verified  (14),  (15),  (17),  (13),  (19),  (20) 
that  both  the  mixing  length  and  the  exchange  coefficient 
vary  over  the  flow  field),  although  certain  modifications 
of  these  ideas  can  produce  reasonably  good  results. 

Prandtl' s  momentum  transfer  theory  is  based  on  the 
hypothesis  that  a  quantity  of  fluid  moves  through  a  stream 
region  as  an  entity,  retaining  its  original  momentum  until 
it  traverses  one  "mixing  length".  At  this  point,  the  fluid 
then  mixes  with  the  surrounding  fluid,  causing  a  local  velocity 
fluctuation  proportional  to  the  difference  between  the 
initial  and  final  mean  flow  velocities.  This  theory  suffers 
from  the  major  drawback  that  it  predicts  similar  rates  of 
transfer  of  momentum  and  energy,  whereas  experimental 
investigations  (14),  (15),  (17),  (18)  show  that  energy  transfer 
is  more  rapid  than  that  of  momentum. 
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The  general  vorticity  transfer  theory  simply 
assumes  vorticity  to  be  conserved  through  the  mixing  process. 
Since  no  practical  calculations  can  be  made  on  this  basis, 

Taylor  proposed  a  modified  vorticity  transfer  theory  which 
is  essentially  the  same  as  the  Prandtl  theory  except  that 
we  are  now  concerned  with  vorticity  instead  of  momentum 
components.  This  theory  does  predict  different  spreading 
characteristics  for  velocity  and  temperature,  but  does  not 
give  prediction  of  profiles  which  are  in  as  good  agreement 
with  experiment  as  does  the  momentum  transfer  theory. 

The  constant-exchange-coefficient  theory  and  its 
modifications  suffer  from  lack  of  physical  interpretation, 
but  seem  to  give  by  far  the  best  agreement  with  experiment  (15). 

In  all  three  theories  the  turbulent  shear  stress 
is  written  T' *  cpjj’—  (by  analogy  to  'T'*  P7lr~  in  the  laminar 
case)  where  P  =  gas  density,  =  mean  velocity  gradient, 

and  -  turbulent  shear  stress.  The  fundamental  question 
concerns  the  nature  of  E  ,  which,  unlike  the  viscosity  v  in 
laminar  flow,  is  not  a  simple  property  of  the  fluid.  Notice 
the  dimensional  equivalence  of  €  to  a  kinematic  viscosity. 

This  quantity  is  described  by  the  theories  in  question  as 
follows : 

-A 

1.  Momentum  transfer  theory:  c  or 

which  makes  r~(fr)x 
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2.  Modified  vorticity  transfer  theory:  <f  -  -.i*  (jr 

Thus  T*  IV  *  However, T  is  still 

proportional  essentially  to  (^.)  ,  with  an  additional 

second  derivative  term  which  is  usually  quite  small. 

pi 

3.  Modified  Exchange  Coefficient  Theory:  £  -  Oc 

where  Uc  is  the  centerline  value  of  the  velocity  at 
a  given  axial  location  and  -  1^(4)  =  the  outer 

momentum  jet  boundary.  The  basic  point  here  is  that 'T' 
is  linearly  proportional  to  .  This  expression, 

simple  as  it  is,  seems  to  provide  best  agreement 

with  experimental  results. 

2.  Analysis 

There  are  three  well-known  analytical  methods 
for  the  study  of  mixing  in  a  turbulent  free  jet.  The  first 
method  is  a  point-source  diffusion  of  momentum,  mass,  and 
temperature.  This  method  is  only  valid  for  large  downstream 
distances  (18)  and  hence  is  not  applicable  in  the  regions  of 
greatest  interest  near  and  around  the  nozzle  exit. 

The  second  method  involves  a  solution  of  the  boundary 
layer  form  of  the  Navier-Stokes  equations,  using  one  of  the 
various  transport  theories  to  relate  the  turbulent  shear  stress 
to  the  other  flow  properties.  These  transport  theories  are: 

(a)  Momentum  transport,  using  the  mixing 
length  concept. 

(b)  Vorticity  transport,  using  a  similar 
mixing  length  concept. 
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(c)  Exchange  coafficient  theory,  using 

a  form  analogous  to  thermal  conduction 
and  diffusion  for  turbulent  momentum 
transfer. 

(d)  Von  Karman  similarity  theory  (21). 

(e)  Statistical  theory  (10),  (21). 

The  third  method  involves  the  solution  of  a  set  of 
integral  equations  of  the  von  Karman  type  (22). 

Of  the  three  methods,  the  last  is  by  far  the  simplest 
and  most  applicable  to  the  present  study.  The  integral 
analysis  leads  to  a  series  of  simultaneous  integro-dif ferential 
equations  which  present  relatively  minor  mathematical 
difficulties  (23),  (24)  compared  to  the  basically  non-linear 
Navier-Stokes  equations.  The  integral  analysis,  of  course, 
suffers  from  the  disadvantage  that  while  differential  equations 
require  only  certain  boundary  values,  the  integral  equations 
require  both  boundary  values  and  specification  of  a  set  of 
initial  profiles  for  the  field  variables.  However,  since  the 
various  profiles  were  measured  with  excellent  accuracy 
(see  data  to  follow)  and  since  integral  analyses  are  generally 
quite  insensitive  to  the  shape  of  the  initial  profiles  (19), 
the  choice  of  the  integral  method  appeared  justified. 

Equations  of  Turbulent  Mixing 

We  now  consider  the  rather  general  problem  of 
the  subsonic  turbulent  mixing  of  a  partially- ionized  gas  with 
a  cool  gas  stream  under  conditions  of  axial  symmetry  as  shown 
in  Figure  1.  The  equations  presented  here  in  final,  non- 
dimensional  form  are  derived  in  Reference  5.  Quantities 
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are  defined  in  the  List  of  Symbols,  Section  V-F 


1.  Argon  -  Argon  Ion  -  Electron  Continuity: 

X  -  I  ^ 


RdR  *7< 


R<iR 


2.  Helium  Diffusion: 
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4 .  Energy: 
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5 .  Spreading  Relations; 
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Assumptions 

The  following  assumptions  were  used  in  the 
analysis.  Each  assumption  is  followed  by  a  statement  of  its 
justification. 

1.  No  Pressure  Gradients:  The  pressure  in  a  free 
jet  of  low  subsonic  flow  is  essentially  constant. 

2.  Equilibrium:  This  assumption  is  necessary  in 
order  to  (a)  calculate  a  temperature  and  (b)  allow  the  use 
of  an  equation  of  state  for  the  multicomponent  gas. 
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Justification  for  use  of  the  equilibrium  condition  is  given 
in  Section  III-C. 

3.  No  External  Body  Forces:  There  were  no  electric 
or  magnetic  fields  applied  to  the  flow,  and  the  effects  of 
gravity  were  certainly  negligible. 

4.  Steady;  Although  not  applicable  on  a  small  time 
scale,  time-averaged  quantities  were  considered  steady.  This 
was  experimentally  verified  by  comparing  integrated 
instantaneous  data  records  with  time-averaged  values. 

5.  Argon  Singly-Ionized  Only;  The  temperatures 
were  much  too  low  for  the  appearance  of  appreciable  doubly- 
ionized  argon  (25),  (26). 

6.  Helium  Does  Not  Ionize;  The  temperatures  were 
also  too  low  for  ionization  of  helium. 

7.  Neutral  Plasma;  Since  the  argon  was  singly- 
ionized  only,  and  there  were  no  fields  present  to  cause 
charge  separation,  the  number  of  positively-charged  ions 
was  equal  to  the  number  of  negatively-charged  electrons  in 
every  region  of  the  flow.  The  Debye  length  was  much  smaller 
than  the  characteristic  dimension  for  temperatures  of  interest. 

8.  Perfect  Gas;  At  the  low  particle  density 
prevailing,  the  accuracy  of  this  assumption  was  at  least 
as  good  as  that  of  the  measurements. 

9.  "Incompressible"  Fluid:  Mach  number  <  0.1. 

Thus  the  test  medium  was  a  variable-density  (due  to  temperature 
and  concentration  changes)  "incompressible"  fluid. 
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10 .  Stagnation  Temperature  Equal  to  Static 
Temperature ;  Again,  Mach  number  <  0.1. 

11.  Prandtl  and  Schmidt  Numbers  Less  Then  Unity: 

This  was  required  by  the  form  of  the  spreading  relations. 
Experimental  data  have  justified  this  assumption.  Also  see 
(14),  (15)  and  Table  3. 

12.  Frozen  Flow  Between  Boundaries:  This  condition 
did  not  actually  exist,  but  because  the  different  boundaries 
were  so  very  close  to  each  other,  as  verified  by  experiment, 
the  effects  of  changes  between  the  boundaries  may  be  neglected. 

13.  Axially  Symmetric  Flow;  This  was  verified 
experimentally. 

14 .  Profiles  in  the  Mixing  Region  were  Similar: 

See  (27)  and  (19).  It  should,  however,  be  pointed  out  here 
that  the  profiles  were  definitely  dissimilar  in  the  so  called 
"potential  core"  near  the  jet  exit,  and  the  similarity 
assumption  was  therefore  not  used  in  the  potential  core. 
Experimental  profiles  in  the  main  mixing  region  were  observed 
to  be  strikingly  similar. 

15.  Inert jal  Frame  of  Reference:  The  boundaries 
which  determine  the  frame  of  reference  were  fixed  (i.e.,  did 
not  accelerate). 

16 .  Profiles  May  be  Represented  by  Cosine  Functions; 
See  (27)  and  (19).  Further,  the  experimentally-observed 
profiles  (see  data)  were  in  reasonable  agreement  with  this 
assumption. 
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17 .  Shear  Work  Term  Negligible  in  Energy  Equation: 

An  order-of -magnitude  analysis  demonstrated  that  this  term 
-3 

was  10  of  the  conduction  terms. 

18.  No  Chemical  Effects  Other  than  Ionization: 

Both  argon  and  helium  are  inert,  monatomic  gases  and  the 
impurity  level  was  too  low  to  require  consideration. 

19.  Radiation  was  neglected  as  a  Mode  of  Energy 
Transfer:  This  is  discussed  in  Section  III-D. 

20.  Reichardt  Theory  was  Used:  As  discussed 
earlier,  this  model  appeared  to  provide  the  best  agreement 
with  experiment  (15). 

As  a  result  of  these  conditions,  the  problem 
essentially  represents  an  extension  of  the  previous  work 
in  the  field  of  turbulent  mixing  in  several  ways: 

(a)  Variations  in  composition,  temperature, 
and  velocity  were  considered  simultaneously. 

(b)  Effects  of  ionization,  very  high  temperature, 
and  extreme  temperature  gradient  were  examined. 

(c)  Plasma-to-gas  heat  transfer  was  significant. 
Nature  of  The  Solution: 

A  detailed  description  of  the  numerical 
solution  on  the  IBM  1620  and  Control  Data  Corporation  1604 
computers  is  given  in  Reference  (5).  The  basic  concepts  of 
this  solution  are  considered  here. 
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Due  to  the  immense  mathematical  difficulties 
associated  with  even  the  simplest  description  of  multi- 
component  turbulent  mixing,  a  closed-form  analytical 
description  is  practically  impossible  without  recourse 
to  rather  unrealistic  assumptions  (23),  (28).  Thus,  a 
method  requiring  a  numerical  solution  on  high-speed  digital 
electronic  computers  was  employed.  While  numerical  solutions 
do  suffer  from  the  difficulty  that  no  direct  functional 
dependence  can  be  obtained,  it  must  be  realized  that  at 
present  the  knowledge  of  turbulent  transfer  processes  at 
very  high  temperature  is  quite  meager  and  realistic  functional 
solutions  do  not  appear  to  be  possible. 

The  fundamental  method  of  solution  of  the  integral 
equations  was  by  successive  Weddle-rule  numerical  integrations 
(29)  of  variables  with  respect  to  the  radial  coordinate  at 
a  fixed  axial  coordinate  (see  Reference  (5)).  Initial 

A 

profiles  (i.e.,  at  the  nozzle  exit,  x  =  0)  must  be  specified  . 
From  this  start,  the  equations  result  in  a  matrix  of  derivatives 
of  the  field  variables  with  respect  to  the  axial  coordinate. 

The  coefficients  of  the  matrix  are  themselves  integrals  in 
the  radial  coordinate  whose  upper  limit  is  one  of  three 
outer  jet  boundaries.  These  boundaries  are  unknown  functions 


These  initial  profiles  were  obtained  from  the  experimental 
measurements.  Other  experimental  inputs  necessary  to  obtain 
the  solutions  appear  in  Table  I. 
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of  axial  position  but  do  obey  the  boundary  condition  that 
they  all  shall  be  unity  (i.e.,  the  non-dimensional  equivalent 
of  the  physical  nozzle  radius)  at  the  nozzle  exit  I  =0. 

With  these  initial  constraints,  the  first  set  of  derivatives 
may  be  calculated,  and  upon  application  of  the  Runge-Kutta 
technique,  new  properties  at  )  =4  are  determined.  There 
are  six  equations  in  the  six  dimensionless  unknowns: 

(a)  Helium  concentration,  CHg 

(b)  Velocity/( Velocity  on  the  centerline  at 

x  =  0),  U 

(c)  Temperature/ (Temperature  on  the  centerline 
and  at  x  =  0 ) ,  © 

(d)  (Outer  concentration  boundary) /(Primary 
nozzle  radius),  R* 

(e)  (Outer  momentum  boundary) /(Primary  nozzle  radius), P 

g 

(f)  (Outer  energy  boundary) /(Primary  nozzle  radius), R0 
Thus  the  solution  is  uniquely  determined  provided  the  equations 
are  independent.  This  solution  at  £  =  then  yields  new 
values  for  the  upper  limits  of  the  integral  coefficients  and 

the  solution  continues  to  propagate  along  the  axial  coordinate, 
continuously  calculating  new  values  for  the  field  variables 
c,  u,  and  ©  .  The  region  of  interest  for  such  solutions  was 
within  R  =  +.  2  and  for  £  (  %  /nozzle  radius)  from  zero  to 
eight.  Beyond  £  =  8  the  jet  was  essentially  decayed. 

The  Matrix 

Upon  differentiating  the  integral  equations 
within  the  integrals  themselves  (using  Leibniz's  Rule)  a 
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resulting  set  of  integro-differential  equations  was  obtained 
of  the  form  (see  Equation  B-47 ,  Reference  (5)). 


av^r)+  (^)+  )♦  o  •  a 


4-o*a}<£ 


where  0  ‘W,(R,R5H($)  Ce-C,(RX)C,C5)  e-eX^O^li) 


Inspection  of  this  matrix  immediately  indicates  a 
non-linearity  due  t">  the  presence  of  the  term  in  (&)(&) 
However,  this  set  is  in  fact  "quasi-linear" ,  that  is,  if  we 
apply  Cramer’s  Rule  to  the  first  three  equations  (i.e.,  striking 
out  the  fourth  row  and  fifth  column)  there  exists  a  set  of 
three  equations  in  four  unknowns.  While  this  is  indeterminate, 
by  using  the  methods  given  in  Reference  (5)  one  may  solve  for 


and  I  in  terms  of  linear  combinations  of 


(-at-) ,  (of-") and  (ar) 

Then,  upon  substitution  of  these  relations  into  the  fourth 
equation,  all  terms  are  linear  i  except  for  the  fourth, 
which  is  quadratic  in  0^^*)  »  Hence(— ^)is  now  the  unknown 
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in  a  simple  quadratic  equation  and  may  readily  be  computed 
in  terms  of  the  various  integral  coefficients  (see  Equation 
B-56,  Reference  (5)). 


Nevertheless 


»■  as(lr> 


frs  now  the  solution  to  an 


algebraic  quadratic  equation,  there  are  in  general  two 
distinct  roots.  However,  as  there  is  dissipative  shear 
between  the  high-velocity  primary  jet  and  the  low-velocity 
secondary  jet,  the  primary  jet  must  accelerate  the  secondary 
flow  at  the  expense  of  its  own  momentum.  Acceleration  of 
the  secondary  flow  to  a  velocity  greater  than  Uc*  ,  which 
was  its  initial  value,  insures  an  outward  spread  of  velocity. 
Since  the  velocity  boundary  is  defined  as  the  locus  of 
points  nearest  to  the  centerline  at  which  U-*  U«.  ,  it  is 
clear  that  in  all  physical  situations  R*  must  increase 
with  axial  position.  Thus  the  only  physically  significant 
solution  is  that  which  has  a  positive  value  for  m-)- 
However,  there  are  still  two  problems  to  be 
considered.  The  first  of  these  is  the  possibility  that  both 
roots  are  positive  numbers.  Fortunately,  this  did  not  occur 
in  any  of  the  numerical  solutions  despite  the  use  of  a  wide 
range  of  input  parameters.  This  does  not,  however,  preclude 
such  a  possibility.  The  physical  interpretation  of  a  doubly- 
positive  set  of  roots  is  not  clear,  and  may  be  worthy  of 
further  investigation. 
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A  second  difficulty  which  did  appear  in  a  few  of 
the  numerical  solutions  concerns  the  discriminant  * 

—  of  the  quadratic  (Equation  B-56,  Reference  (5)). 

If  the  root  is  imaginary,  and  )  »  and  therefore 

all  the  other  axial  derivatives  of  the  field  variables, 
become  complex  numbers  with  an  imaginary  part.  Since  such 
results  (i.e.,  complex  values  of  velocity,  temperature  and 
concentration)  appear  to  have  no  direct  physical  significance, 
such  solutions  were  considered  meaningless.  The  physical 
interpretation  of  these  solutions  may  also  bear  further  study. 

The  stability,  or  presence  of  "meaningful"  solutions 
in  the  sense  of  the  above  discussion,  is  related  to  the 
integral  coefficients  which  are  correspondingly  related  to 
various  input  parameters.  Since  the  purpose  of  the 
analytical  study  was  primarily  to  furnish  information  on 
the  nature  of  the  turbulent  mixing  process,  a  complete 
investigation  of  "meaningless"  or  "unstable"  solutions  has 
not  been  performed.  However,  certain  broad  and  rather 
significant  conclusions  may  be  drawn: 

1.  The  boundaries  must  be  reasonably  close  to 
each  other.  Whenever  large  values  of  A^  or  AjWere  used 
(corresponding  to  boundaries  greatly  separated  at  reasonable 
axial  positions)  solutions  were  of  the  "unstable"  type. 

This  might  well  be  the  result  of  severe  violation  of  the 
assumption  of  frozen  flow  between  the  boundaries. 
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2.  Large  values  of  the  forcing  functions;  i.e., 
small  values  of  the  equivalent  turbulent  Reynolds  numbers, 
resulted  in  "unstable"  solutions. 

3.  Excessive  dominance  of  any  one  transfer 
process  (i.e.,  material  diffusion,  momentum  transport,  or 
turbulent  heat  transfer)  expressed  via  the  driving  parameter, 
resulted  in  "unstable"  solutions.  Specifically,  "unstable" 
solutions  were  found  to  occur  whenever 

(a)  ft,  ■>  O.  0080 
0,0100 

<b>Nfeo»  W*eNp*=  Kl*  -200 

^  <  *0£> 

^Re  Np«  =  C  200 

(c)  Cases  in  which  any  one  of  the  parameters 

differed  by 

more  than  a  factor  of  10  from  any  other. 

The  accuracy  of  a  numerical  solution,  particularly 
those  involving  repeated  numerical  integration  schemes  and 
the  Runge-Kutta  method, is  always  subject  to  question.  The 
best  method  of  checking  the  extent  or  at  least  the  nature 
of  the  error  is  by  using  finer  and  finer  increments.  It  was 
found  that  variation  of  increment  or  "mesh"  size  over  an 
order  of  magnitude  had  no  effect  on  the  solution  beyond  the 
fifth  digit  after  the  decimal  point.  This  excellent  accuracy 
was  probably  due  to  the  smooth,  monotonic,  single-valued 
nature  of  all  the  various  functions  employed. 
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Finally,  the  assumption  of  "similar  solutions" 

(i.e.,  functionally  similar,  affinely  related  profiles) 
for  the  field  variables  must  be  considered.  The  work  of 
Squire  and  Trouncer  (27)  indicates  that  the  solutions  are 
similar  in  the  mixing  region,  but  are  not  similar  in  the 
potential  core.  This  was  observed  experimentally.  Thus, 
in  order  to  obtain  a  more  physically  reasonable  solution 
it  was  necessary  to  include  a  potential  core.  However,  as 
the  number  of  available  equations  had  been  exhausted, 
it  was  necessary  to  do  this  empirically.  Reduced  experimental 
data  provided  the  necessary  information.  This  somewhat 
arbitrary  constraint  was  then  included  in  the  program  for 
the  numerical  solution.  Comparison  of  results  with  and 
without  a  potential  core  indicated  that  there  was  little 
or  no  effect  on  spreading,  although  there  was  observed  the 
anticipated  effect  upon  both  initial  radial  profiles  and 
initial  axial  decay.  To  provide  a  bit  more  generality, 
the  length  of  the  cone-shaped  core  was  allowed  to  vary. 

This  variation  had  little  effect  on  the  nature  of  the 
radial  spreading  of  velocity,  concentration,  and  temperature. 

The  numerical  results  of  this  analysis  and  their 
comparison  with  experimental  data  are  given  later  in 
Section  IV. 

B.  Laminar  Transport  Properties 

Considerable  theoretical  study  and  analysis  of 
the  transport  properties  for  helium  and  ionized  argon 
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mixtures  was  performed  during  the  contract  period.  This 
work  is  shortly  to  be  completed,  and  a  technical  report 
presenting  the  analysis  and  results  will  be  issued  at 
that  time.  An  introductory  statement  of  the  problem  and 
the  methods  of  attack  are  presented  here.  Figures  2  and  3 
show  the  results  of  viscosity  calculations  (typical  of  the 
general  class  of  transport  property  results  to  appear  in  the 
forthcoming  technical  report.) 

1.  The  Problem 

The  problem  considered  here  is  the  determination 
of  the  viscosity,  thermal  conductivity,  diffusion,  and 
thermal  diffusion  coefficients  of  argon-helium  mixtures 
at  temperatures  up  to  15,000°K  and  pressures  around  one 
atmosphere.  At  the  high  temperature  end  of  this  region,  a 
considerable  fraction  of  the  argon  is  singly  ionized.  The 
helium,  having  a  much  higher  ionization  potential,  can  be 
considered  as  inert.  The  fraction  of  doubly  ionized  argon 
is  negligible. 

2.  Composition 

Since  the  mole  fraction  of  argon  atoms  or  ions 
with  excited  electronic  states  is  always  small  compared  to  the 
mole  fraction  in  their  ground  states,  the  composition  of 
the  mixture  could  be  considered  as  a  four  component  system 
consisting  of  the  following: 
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(a) 

Argon  atoms 

(ground  state) 

* 

(b) 

Argon  ions 

(singly  ionized  in  their 

ground  state) 

i 

(c) 

Electrons 

(d) 

Helium  atoms  (ground  state) 

The  system  was  also  considered  to  be  in  thermal 
equilibrium,  as  will  be  discussed  in  III-C.  The  mole 
fractions  of  argon  neutrals,  argon  ions  and  electrons  was 
obtained  from  the  statistical  mechanics  computation  of  (30) 
and  (31)  using  the  expressions 

^  =  2.309  x  10~6 |lj  T2  *  5  51exp 

where  P  =  pressure  in  atmospheres 
T  =  temperature  in  °K 

Xg  =  mole  fraction  of  electrons  or  argon  ions 
=  mole  fraction  of  argon  neutrals 
This  is  essentially  the  Saha  equation  with  empirical 
corrections  for  the  effect  of  the  electronic  partition  functions. 
It  should  be  pointed  out  that  an  uncertainty  in  the  results 
exists  because  of  the  arbitrary  position  of  the  cutoff  in 
the  argon  electronic  partition  function.  This  uncertainty  is 
not  large,  however,  and  is  not  felt  to  appreciably  affect  the 
subsequent  calculation  of  the  transport  properties. 

3.  Transport  Property  Theories 
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The  difficulties  encountered  in  the  calculation 
of  transport  properties  can  be  grouped  into  three  types: 

(a)  The  obtaining  of  a  suitable  theory. 

(b)  The  obtaining  of  the  physical  information 
needed  by  the  theory;  e.g.,  collision  cross  sections. 

(c)  Carrying  out  the  complex  calculations. 

The  theories  for  the  computation  of  transport 

properties  can  be  divided  into  two  groups,  the  formal  or 
"exact"  theories  and  the  admittedly  approximate  ones.  For 
any  exact  theory,  and  probably  for  the  approximate  ones,  the 
labor  of  computing  the  transport  properties  of  multicomponent 
mixtures  is  so  enormous  that  the  use  of  electronic  computers 
is  essential.  All  the  computations  were  performed  on  an  IBM 
1620  digital  computer. 

Obtaining  the  physical  information  used  by  the 
theory  is  usually  the  most  difficult  part  of  the  entire 
procedure.  Since  the  information  is  usually  fragmentary 
or  unavailable,  considerable  personal  judgment  enters  into 
this  phase.  As  an  illustration  of  the  possible  difficulties 
in  the  interaction  of  an  argon  atom  with  an  argon  ion  at 
energies  corresponding  to  the  temperatures  of  interest 
(1  to  2  electron  volts),  the  long-range  polarization  force 
is  not  important,  but  it  is  nevertheless  considered  so  in 
certain  references. 
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In  the  computation  of  the  transport  properties 
of  an  ionized  or  partially  ionized  gas  the  difficulty  of 
obtaining  a  suitable  theory  arises.  For  neutral  rare 
gases  at  temperatures  and  pressures  where  they  are  nearly 
perfect  gases,  the  transport  properties  can  be  rather 
accurately  computed  by  the  use  of  the  Chapman-Enskog  theory 
(32,  33).  The  mechanics  of  using  this  theory  have  been 
so  developed  (33)  that  its  use  by  engineers  in  accurate 
computation  of  transport  properties  is  nearly  universal. 
Unfortunately  in  applying  this  technique  to  ionized  and 
partially  ionized  gases,  problems  arise  which  at  the  very 
least  make  the  computations  more  difficult,  and  perhaps  throw 
the  reliability  of  the  entire  method  in  doubt. 

The  Chapman-Enskog  method  is  derived  from  the 
Boltzmann  equation  in  which  the  collision  term  considers 
only  binary  collisions.  It  is  assumed  that  during  most 
of  the  time  a  particle  exists  it  is  far  enough  from  the 
other  particles  that  the  effect  of  their  intermolecular 
force  fields  on  the  trajectory  of  the  particle  is  negligible. 
Then  events  occur  in  which  two  particles  pass  so  close 
together  that  we  can  say  an  encounter  or  collision  has 
occurred.  It  is  assumed  that  the  likelihood  that  three  or 
more  particles  will  be  close  enough  to  significantly 
influence  each  other's  paths  is  negligible  compared  to  the 
two-body  collisions.  The  derivation  of  the  Chapman-Enskog 
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theory  leads  to  the  "collision  cross  sections"  and  "collision 
integrals"  by  which  the  specific  intermolecular  forces  enter 
into  the  computations;  When  the  coulombic  potential  is 
introduced,  the  collision  cross  sections  and  collision 
integrals  become  divergent  as  the  integration  is  extended  to 
large  "impact  parameters"  (the  minimum  distance  the  two 
particles  would  approach  one  another  if  their  paths  had 
not  been  deflected).  This  is  overcome  in  (32)  by  cutting 
off  the  integration  at  the  mean  interparticle  distance. 

It  is  probably  more  accurate  to  use  the  Debye  length  as 
the  cutoff. 

This  difficulty  is,  however,  symptomatic  of  the 
unusual  properties  of  the  Coulomb  potential.  This  potential 
falls  c  ~f  so  slowly  compared  to  the  usual  intermolecular 
potential'  that  the  cumulative  effects  of  long-range  weak 

S 

interactions  on  the  trajectory  of  a  charged  particle  are 
greater  than  the  effect  of  "collisions"  (34,  35).  The 
differences  in  the  behavior  of  plasmas,  even  without 
externally  applied  electric  and  magnetic  fields,  from  that 
of  ordinary  gases  is  so  great  that  a  new  kinetic  theory  is 
really  needed  to  deal  with  them.  Since  this  has  not  been 
developed  to  the  point  that  it  can  be  used  to  compute 
transport  properties,  ordinary  kinetic  theory  with  suitable 
alterations  was  used. 
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A  second  difficulty  with  the  Chapman-Enskog 
method  is  that  almost  all  computations  use  only  the  lowest 
approximation  in  the  variational  scheme  of  (33).  This  is 
usually  sufficient  to  give  excellent  results,  but  for 
plasmas,  at  least  a  second  approximation  should  be  made  for 
the  diffusion  coefficients  (36,  37). 

In  an  ordinary  gas  the  only  important  characteristic 
distance  is  the  mean  interparticle  distance.  In  a  plasma 
the  important  characteristic  distances  are 

(a)  The  mean  impact  parameter  for  a  90°def lection,! 

(b)  The  mean  interparticle  distance  d. 

(c)  The  Debye  length  X  . 

For  most  plasmas  of  interest, 

p  <<  d  <<  X 

*o 

Out  to  pQ  the  collisions  are  of  the  strong-encounter,  two- 
body  type  considered  in  ordinary  gases.  Beyond  the  Debye 
length  the  effect  of  a  particle  is  assumed  to  be  shielded, 
and  particles  beyond  are  assumed  not  to  have  an  effect  in 
the  computation  of  transport  properties  (34).  The  important 
regions  are  between  the  Debye  length  and  the  mean  interparticle 
distance  and  between  this  distance  and  pQ.  To  handle  this 
weak  interaction  successfully,  a  model  has  been  developed  which 
uses  instead  of  the  ordinary  Boltzmann  equation  with  a  binary- 
encounter  collision  term,  a  Boltzmann  equation  with  a  special 
collision  term  supposedly  more  suitable,  called  the  Fokker- 
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Planck  equation  (34,  35).  This  has  been  used  thus  far  to 
compute  the  transport  properties  of  fully  ionized  gases. 

When  the  results  are  compared  to  those  of  the  Chapman- 
Enskog  theory,  suitably  carried  out  to  enough  approximations 
(37),  similar  results  are  found. 

There  are  numerous  approximate  schemes  for 
computing  transport  properties  based  on  mean  free  path 
concepts,  etc.,  which  may  or  may  not  use  refinements  like 
persistance-of-motion  corrections.  It  is  generally  known 
that  these  methods  can  be  fairly  successful  for  viscosity 
and  thermal  conductivity  calculations,  but  are  poor  for 
diffusion.  A  "Relaxation  Time"  method  was  used  in  the 
present  analysis  for  the  computation  of  viscosity.  The 
results  agree  fairly  well  with  results  of  the  complex 
Chapman-Enskog  method,  as  can  be  seen  from  Figures  2  and  3. 
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C.  Arc jet  Equi I ibr ium 
I  .  Introduction 

The  state  of  a  fluid  in  a  region  can  not  be  in  equilibrium 
if  there  are  finite  gradients  of  temperature,  composition,  or  velocity. 

The  equilibrium  state  for  a  fluid  in  a  small  region,  nonetheless  containing 
many  molecules,  is  that  state  the  fluid  would  eventually  approach  were 
it  suddenly  isolated.  For  a  great  number  of  problems  the  actual  values 
of  quantities  such  as  density,  mole  fractions,  translational  energy. 
Internal  energy,  etc.  are  negligibly  different  from  equilibrium  values 
so  that  these  flows  are  called  "equilibrium"  flows.  One  could  relax 
the  criteria  to  ell  flows  in  which  the  differences  from  equilibrium 
were  at  most  only  a  few  percent  from  equilibrium  values.  The  latter 
criterion  is  the  one  used  in  this  work. 

2.  Problem 

The  problem  considered  here  is  whether  the  ergon  plasma 
generated  in  o  plasma jet  under  certain  conditions  is  in  equilibrium 
as  it  leaves  the  nozzle  of  the  plasmajet,  and  if  it  remains  in 
equilibrium  as  it  is  rapidly  cooled  in  the  mixing  region  downstream. 

The  flow  is  at  low  subsonic  Mach  numbers,  and  thus  is  at  ambient  (in 
this  case  atmospheric)  pressure.  The  peak  temperature  of  the  jet 
leaving  the  nozzle  is  about  I4,000°K. 

Two  modes  of  operation  are  considered,  turbulent  and 
laminar.  In  the  former  case  peak  flow  velocities  are  about  500  ft/sec, 
mass  flows  around  4  x  I0“5  lb/ sec,  current  about  2000  amp,  and  voltages 
about  30  volts.  The  nozzles  used  had  a  length  of  eight  Inches  from 
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the  ere  to  the  nozzle  exit.  In  the  laminar  case  the  velocities 
ere  much  lower*  the  flow  rates  about  5  x  10”*  lb/ sec*  the  current 
ebout  500  amp,  end  the  voltages  about  25  volts.  The  nozzles  used 
are  about  I  1/2  inches  long.  Thus  the  times  spent  in  traveling 
between  the  arc  and  the  nozzle  exit  are  about  the  same  in  both 
cases.  The  gradients  in  the  mixing  region  will  be  more  severe  in 
the  turbulent  flow;  hence  "equl  Hbrium"  of  the  turbulent  Jet  is 
the  limiting  case.  Equilibrium  of  the  turbulent  Jet  was  investi¬ 
gated  in  Ref.  (2).  This  analysis  is  an  Improvement  of  the  treat¬ 
ment  used  therein. 

For  this  problem  the  criteria  for  equilibrium  can 
be  stated (39)  as: 

Are  the  mole  fractions  of  argon  neutrals,  argon  ions, 
and  electrons  near  equilibrium? 

Are  the  electron  and  Ion  temperatures  close? 

Are  the  populations  of  the  excited  electronic  states 
near  equilibrium? 

The  concept  of  different  electron  and  heavy-pert Ic I e 
temperatures  results  from  the  extreme  ineffectiveness  of  energy 
transfer  In  electron-heavy  particle  elastic  collisions  compered 
to  that  in  collisions  between  particles  of  equal  mess.  The  electrons 
rapidly  attain  a  near-tiaxwel I  Ian  distribution  after  a  few  electron- 
electron  collisions,  and  the  heavy  particles  also  attain  a  near- 
Maxwell  ian  distribution  after  a  few  heavy  particle-heavy  particle 
collisions.  However,  the  two  distributions  can  correspond  to  different 
temperatures,  called  the  electron  and  the  heavy-particle  temperatures 
respectively.  It  takes  the  order  of  40,000  electron-heavy  particle 
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collisions  to  "thermal ize"  an  electron. 

3.  The  Physical  Processes 

Energy  Is  added  to  the  gas  as  it  passes  through  a 
DC  arc  struck  between  a  tungsten  cathode  and  a  copper  anode.  Both 
electrodes  ere  internally  water  cooled.  Almost  the  entire  current 
consists  of  electrons.  It  is  believed  that  the  electrons  are 
thermion  leal ly  emitted  from  the  cathode.  Note  that  a  current  of 
2000  amperes  corresponds  to  about  10^  electrons/ sec. 

The  energy  of  the  arc  first  goes  into  Increasing 
the  electron  fra ns  I at iona I  energy.  Many  ions  may  be  formed  by  the 
electrons  in  the  arc  and  many  are  formed  by  thermal  ionization 
downstream  of  the  arc.  At  argon  flow  rates  of  4  x  I0”5  lb/ sec 
the  number  of  ions  that  must  be  created  tor  a  I5,000ft<  jet  In 
equilibrium  is  2.6  x  IQ^/sec.  It  Is  clear  that  because  there  Is 
alwoys  such  a  large  number  of  electrons  present,  atom-atom  ionization 
will  be  negligible  compared  to  electron-atom  Ionization.  It  will  be 
shown  that  thermal  ionization  by  electrons  Is  sufficient  to 
produce  the  required  number  of  ions. 

The  cross  section  for  the  ionization  of  a  ground  state 
argon  atom  by  an  electron  Is  fairly  well  known.  Fox  (40)  and  Tozer 
o  ,d  Craggs  (4i)  show  a  roughly  linear  rise  In  cross-section  from  a 
threshold  at  the  ionization  potential.  Petschek  and  Byron  ( 39)  use 
a  cross-section  with  a  threshold  at  the  critical  potential,  11.5 
electron-volts.  This  they  attribute  to  a  multiple-stage  Ionization 
process.  The  lower  number  of  atoms  in  the  excited  states  Is  balanced 
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by  the  increased  number  of  electrons  capable  of  ionizing  the  excited 
states  as  compared  to  the  ground  state.  The  cross  section  for 
ionization  of  an  excited  state  is  probably  higher  than  that  for  the  ground 
state  (42).  For  an  optically  thin  plasma  the  excited  atom  populations 
could  be  much  smaller  than  the  equilibrium  values  because  of  the 
relative  decrease  In  radiative  excitations  as  compared  with  radiative 
decay  of  the  excited  states.  In  this  case,  colllslonal  excitation 
must  be  sufficient  so  that  the  spontaneous  radiative  decay  does  not 
depopulate  the  excited  levels.  Typical  lifetimes  for  the  excited 
states  before  radiative  decay  are  10"®  to  10"®  sec.  (43).  This  Is 
short.  However,  there  are  at  least  two  mechanisms  which  increase  the 
effective  lifetimes:  first,  orgon  has  meTastable  excited  states,  and 
second,  under  the  conditions  of  the  problem  the  "mean  free  path" 
of  resonance  radiation  is  extremely  small.  Thus  the  energy  must  be 
passed  from  atom  to  atom  before  leaving  the  plasma,  and  it  therefore 
appears  likely  that  most  ionization  will  occur  via  at  least  a  two- 
step  process. 

The  ionization  process  removes  15.75  electron-volts 
per  ion  from  the  translational  energy  of  the  electrons.  It  is 
reasonable  to  believe  that  the  bulk  of  the  electrons  will  not  have 
energies  above  the  critical  potential  of  argon,  11.5  electron-volts, 
due  to  Inelastic  collisions,  and  hence  it  is  quite  conservative  to 
consider  that  electrons  lose  their  translational  energy  (from  11.5 
electron-volts  to  thermal  energies)  only  by  elastic  collisions. 

The  electron-ion  cross  section  is  so  much  larger  than  the  electron- 
atom  cross  section  that  we  may  ignore  energy  transfer  by  the  latter- 
type  collisions.  It  has  been  ccmputed(44)  that  at  only  0.1# 
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Ionization,  ol^ctron-ion  energy  transfer  by  elastic  col  Melon  Is 
equal  to  the  electron-atom  energy  transfer. 

In  the  mixing  region  of  the  case  under  consideration, 
energy  Is  lost  by  the  heavy  particles  to  the  atoms  of  a  cool  outer 
helium  flow.  The  electrons  must  transfer  their  energy  to  the  atoms 
or  ions  in  order  to  keep  the  electron  temperature  close  to  the 
heavy  particle  temperature.  As  the  temperature  falls,  the  equilibrium 
concentration  of  electrons  falls,  requiring  that  the  electrons 
recombine  with  the  ions.  Unfortunately,  even  today  the  subject  of 
electronic  recombination  in  rare  gases  Is  not  well  understood. 
Physicists  have  devoted  a  great  deal  of  effort  to  this  subject, 
but  all  too  often  their  results  confuse  rather  than  clarify. 

As  an  illustration  of  the  difficulties,  20  years  ago  the  accepted 
rate  of  recombination  in  argon  was  1/1000  the  presently  accepted 
value.  The  rate  of  recombination  is  usually  given  by 

=  _  =  ,  quasi  neutral  I  ty  condition). 

*Lt 

The  recombination  coefficient,  «<  ,  Is  a  strong  function  of  the 

temperature  and  also  a  function  of  the  electron  density.  Most  of 
the  measurements  of  at  were  at  low  temperatures  and  low  electron 
densities.  Drastic  extrapolations  ore  needed  to  obtoin  values 
for  plasmajet  conditions. 

Even  if  accurate  experimental  values  of  the 
recombination  coefficient  for  the  conditions  of  interest  existed. 

It  would  still  be  of  great  value  to  know  the  type  of  mechanisms 
predominating.  If  the  energy  released  by  the  recombination  does 


42 


not  go  back  into  alactron  kinetic  energy,  but  is  either  lost  to 
the  plasma  by  radiant  emission  or  goes  into  heavy  particle  kinetic 
energy,  the  electron  temperature  will  more  closely  follow  the 
heavy  particle  temperature. 

The  presently  accepted  recombination  mechanisms  for 
argon  aret  (a)  radiative  recombination,  (b)  dielectronlc  recombination, 
(c)  dissociative  recombination,  end  (d)  three  body  recombination. 
Radiative  recombination  Involves  the  collision  of  an  electron  and  a 
positive  ion  with  the  radiative  emission  of  energy.  This  mechanism 
has  been  studied  using  quantum  mechanics.  For  the  hydrogen  atom, 
tables  of  recombination  coefficients  exist,  together  with  a  break** 
down  of  the  relative  importance  of  the  recombination  to  the  various 
energy  levels  of  the  atom,  (45,  46).  They  show  that  in  hydrogen 
large  fractions  of  the  recombination  are  to  the  ground  state  and 
the  other  low  levels.  Values  oV  The  radiative  recombination 

rate  «4r*  appear  below 


RADIATIVE  RECOMBINATION  COEFFICIENT  «r  IN  HYDROGEN  (8,  9) 
alMDLiS&I  Urirm  3/ 


8,000 

16,000 

32,000 


4.83  x  10-13 
2.93  x  I0’15 
1.72  x  IO-«5 


From  the  results  given  In  (45)  and  (46),  it  appears  the  recombination 

coefficients  for  the  rare  gases  are  close  to  that  for  hydrogen.  If 
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this  Is  the  case,  then  radiative  recombination  will  be  a  small  but 
not  negligible  part  of  the  total  recombination. 

□(electronic  recombination  involves  the  collision  of  an 
electron  and  a  positive  ion  in  which  the  two  ere  bound  in  a  level 
above  the  ionization  potential.  This  may  occur,  for  example, 
by  excitation  of  two  electrons.  In  order  for  the  recombination 
to  proceed,  this  level  must  be  deexclted  by  either  radiative 
emission  or  by  superelastic  collisions  before  It  reionizes. 
Quantitative  models  for  di electronic  recombination  are  very  crude, 
but  it  appears  that  the  rates  efe  at  most  as  large  as  those  of 
radiative  recombination.  At  high  temperatures  and  high  electron 
densities  dielectronlc  recombination  Is  probably  negligible. 

In  the  past  fifteen  years  experiments  in  which  cold 
argon  at  low  pressures  was  ionized  by  microwave  pulses  semeed  to 
indicate  that  the  dominating  mechanism  was  dissociative  recombination 
(52)-  The  existence  of  rore  gas  molecular  Ions  has  been  verified. 
According  to  this  model,  the  recombining  electron  dissociates 
the  molecule,  transferring  its  energy  into  breaking  the  molecular 
bond,  translation  of  the  atoms,  and  electronic  excitation.  Most 
of  the  energy  probably  appears  os  electronic  excitation. 

In  order  that  this  mechanism  be  the  dominant  mechanism 
for  recombination  in  a  near-equl 1 ibrium  plasma,  the  inverse  process, 
associative  ionization, 

A  *  A  — *  At  +  e 

must  be  the  dominant  mechanism  of  ionization,  in  accordance  with 
the  principle  of  detailed  balancing.  This  mechanism  has  not  been 
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reported  to  bo  Important  In  raro  gesos.  If  mo  const  dor  that  tho 
formation  of  A £  ,  such  as  by  three  body  collisions  of  the  type 

A t  +M 

Mhoro  M  is  soma  third  body* 

ere  vary  rapid  mo  can  writo  an  equilibrium  rotation  for  < 

,  - 

,T  Na  Kt*  k 

To  compute  tho  value  of  tho  equilibria  constant,  wo  need  tho 

value  of  tho  binding  energy  end  othar  molecular  constants.  Tho 

binding  onergy  for  Htt  is  about  2. Ir.y. Ctl),  so  If  wo  assume  that 

tho  formation  of  Ax  is  comparable,  wo  take  Its  binding  energy  at 

about  2  e.V..  It  can  be  shown  that  a  I  though  i  s  very  small  at 

high  temperatures.  It  con  bo  quite  largo  at  low  tmaporaturos. 

Hence  in  tho  microwave  experiments,  since  K  ^and  He*  (compared 

to  Ha  )  are  large,  Na£  can  be  large,  and  dissociative 

recombination  will  dominate.  At  plasmajet  conditions,  however, 

the  mole  fraction  of  Ax  is  estimated  at  10*^  to  I0-6, 

Although  smell,  this  mole  fraction  is  not  so  small  that  dissociative 

recombination  can  definitely  be  neglected.  If  the  cross  section 

for  dissociative  recombination  of  At  Is  large,  this  mechanism 

can  still  be  important.  The  uncertainties  are  thus  so  greet 

that  it  can  not  be  determined  if  dissociative  recombination  in 

argon  at  plasmajet  conditions  is  negligible,  significant,  or 

dominant. 

In  the  past  few  years  there  has  been  a  great  deal  of 
work  on  the  study  of  three  body  recombination 

e+e+A*  ?=*  A+# 
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For  th*  hydrogen  atom  (46)  gives  the  recombination  coefficients 
listed  In  below. 


THREE  BOOY  RECOMBINATION  COEFFICIENT  FOR  HYDROGEN 


(cm.5/ sec.) 


ELECTRON  DENSITY  (cm.-5) 

TEMPERATURE  <°K> 

8000 

16,000 

32,000 

I0»6 

,  8.4  x  I0*U 

5.0  x  I0*12 

7.3  x  I0“13 

.o'7 

3.4  x  I0‘,& 

1.4  x  I0*M 

1.8  x  I0-12 

I018 

2.5  x  I0”9 

9.6  x  I0"M 

1.2  x  I0"M 

Extrapolations  of  value*  for  heliuai  (47)  give  PC  %  I0~' ?  cm.^/sec 
for  He  •  I017  ,  T  ■  15,000  K,  which  is  about  the  same  as  that  for 
hydrogen.  The  three-body  recombination  coefficient  tor  argon  is  also 
probably  close  to  the  values  quoted  in  Table  II  for  hydrogen.  Existing 
models  indicate  that  recombination  directly  to  the  ground  or  low-lying 
states  Is  rare,  but  that  most  of  th*  recombinations  take  place  to 
intermediate  levels.  The  inverse  of  this  process,  electron-atom 
collision*!  ionization,  is  generally  accepted  to  be  th*  dominant 
ionization  process.  It  Is  most  likely  thet  In  th*  present  region 
of  interest,  this  recombination  mechanism  is  th*  dominant  on*. 

By  standards  used  in  recombination  studies,  th*  jet 
exhaust  must  be  considered  as  full  of  Impurities.  Although  th* 
rare  gases  have  the  lowest  recombination  coefficients,  impurities 
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Mill  accelerate  the  recombination.  At  low  temperatures  (5000°K)  the 
Impurities  act  as  "seeding"  materials  tor  the  production  of  electrons, 
and  thus  the  electron  concentration  could  be  orders  of  magnitude 
higher  than  the  equilibrium  value  tor  argon.  However,  the  mole 
fraction  of  the  electrons  appears  to  be  sufficiently  small  that 
their  effect  on  the  thermodynamic  properties  will  not  be  significant. 


4.  Elastic  Col  I isions 


The  problem  of  the  rate  of  energy  transfer  by  elastic 
collisions  would  be  straightforward  were  it  not  for  the  unusual 
properties  of  the  Coulomb  potential.  Because  of  the  slow  felloff 
of  this  potential,  the  effects  of  small  deflection,  long  range 
encounters  are  more  Important  than  those  of  close  encounters. 

Energy  transfer  by  electron-atom  collisions  is  negligible  compared 
to  electron-ion  Interactions  if  there  is  even  \%  ionization. 

Since  the  Important  collisions  are  no  longer  of  the  two-body  type, 
the  Fokker-Plenck  equation  can  be  used  to  treat  this  case.  The 
results  are  essentially  the  same  as  those  of  ffctschek  end  Byron  (39). 

fVtschek  and  Byron's  method  consists  of  taking  an 
expression  tor  the  energy  transfer  occurring  in  an  elastic  electron-ion 
collision  of  a  given  deflection,  and  suitably  Integrating  over  all 
coil isions 439) .  The  electrons  are  assuaied  to  be  in  a  Maxwellian 
distribution,  and  the  electron  mess  is  considered  small  coaipere« 
to  the  ion  mass.  The  result  (37)  is 


At  *  tfiiii  e* 

Mi 
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where  energy  transfer  per  unit  uolUM 

fje  ■  electron  particle  density 
N;  •  ion  particle  density 
Ms  ■  electron  mass 
M*  •  Ion  mass 
Jl  *  Boltzmann  constant 
■  electron  temperature 
Tt  •  ion  temperature 
®  ■  electron  charge 

(C.  G.  S.  electrostatic  units  ere  used  throughput  this  report).  The 
above  formula  can  be  expressed  by 

cHe*(  lkl±) 

Te  '  . 

where  C  Is  almost  constant  at  5  x  10"**  erg,  os.'/nc. 

Turning  now  to  the  method  using  thd  Fokkor -Planck  equation, 
assuming  special  homogeneity  and  no  body  forces,  the  Boltzmann 
equation  for  the  electrons  is 

9fr  _  electron-electron  collision  tons  t  eloctron-lon 
gt  collision  term  +  electron-otpo  solllslon  term 

Assuming  a  Maxwellian  distribution  of  electrons*  the  electron-electron 

collision  term  is  zero.  We  neglect  the  slectren-eten  collision  term 

for  even  moderate  degrees  of  ionization  as  discussed  per  Her.  Using 


a  Fokker -Planck  form  for  the  electron-ion  collision  tens 

TT  "N  4  f*  V‘  A  '  M  ‘ F'  Y W  j  -  Mr) v, 

where:  1  •  unit  tensor  —  %Jti  ’  •  _ 


<f  ■  +7Tc*  rtu4.fr  cl*  JL  A.  •  ertss  section 

at**  He1  M/* 


2  fa+N, 
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A  -- 


C’f+rrve 


Vr 


relative  velocity 


Multiplying  by  1/2  Mf  tTtf1  eft#  Integrating 
over  all  electron  velocities  «e  obtain  after  considerable 
mathematical  effort 


This  result  is  essentially  the  saaw  ee  that  of  (39X  The  alight 
difference  In  the  logarithm  la  due  to  a  difference  In  definition 
of  the  Debye  length. 

The  enthalpy  per  unit  volume  at  the  electrons  can  b# 
..pc..,**,.  w#.jrN,iT#+FiN, 

eel  9 

where  £L  *  ■  Ionization  energy  per  electron  •  25.2  a  10  erg 


ut 


J.*  zrrzr***1* 

We  obte’-  U,  He  _ 


* 


s[i  Ne  J  ^ r«  ♦  M’eJ 


The  first  term  represents  the  leas  In  thermal  energy  aauaed  by  the 
cooling,  the  second  term  repreaenta  the  lead  In  thermal  energy 


caused  by  the  removal  of  seme  of  the  free  electrons  by  reeeami nation 
and  the  third  term  represents  the  leas  ef  Inn  I  let  lea  energy  e#  these 
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electrons.  The  relative  Magnitude  of  that*  farms  Is  roughly  In 
tha  ratio  i':4il6.  From  tha  experimental  data  of  (5),  estimates 
of  ‘jtf?  can  be  made;  s.g., 

•  velocity  x  gradient  of  tsmperaturo  along  a  streamline 
Taking  maxiMum  values  from  (15),  we  select  a  speed  of  15,000  cm. /sec. 
(500  ft. /sec.)  and  a  temperature  gradient  of  l5,OOOPK/cm.  This 
results  In  a  value  tor  of  2.3  x  10®  °K/sec.  which,  tor 

lack  of  more  precise  information,  will  be  used  as  an  extreme 
upper  bound.  The  r«*te  of  electron  en+halpy  decay  Is  given  in 
below  as  a  function  'f  the  rate  ot  Temperature  decay  and  local 

electron  temperature  (the  numbers  in  parentheses are  the  individual 
contribution  ot  the  three  farms).  Note  that  I  kw/cm3  ■  10*®  erg/cm5  * 


RATE  OF  ELECTION  ENTHALPY  DECAY 


5  x  I07 
I  x  10® 
2.3  x  10® 


T  •  I2,00CPK 


Enthalpy  Decay  (kw/cm3) 


T  •  I3,00tf*  T  •  U,00<fl< 


(.02  ♦  .17  ♦  1.02)  (.04  ♦  .19  ♦  1.02) 


12.10 

27. SO 


12.40 

28.30 


12.00 

29.00 


The  problem  is  now  to  relate  the  electron  enthalpy  decay 
rate  a?  to  Af  •  For  example.  If  there  were  no  net  energy 
transfer  by  Inelastic  collisions,  the  two  would  be  equai,  whereas 
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if  the  recombination  process  resulted  In  the  energy  of  the 
electron  going  entirely  into  heavy  particle  translational  energy, 
thenAEwould  he  roughly  equal  to  the  first  term  in  the  expression 
♦or  fL.P* 

Jit, 

5.  Energy  Exchange  Processes  in  a  Plasma 

To  find  the  fraction  of  the  enthalpy  decay  of  the  electrons 
which  mtsrt  be  transferred  by  elastic  collisions,  a  qualitative 
understanding  of  the  processes  that  take  place  In' the  plasma  as  It 
cools  Is  necessary,  as  well  as  some  quantitative  estimates  of  the 
relative  Importance  of  the  various  processes. 

Although  our  knowledge  of  these  processes  is  very  incomplete, 
and  It  is  therefore  not  certain  that  any  model  we  may  postulate  will 
be  realistic,  it  Is  nevertheless  of  interest  to  at  least  speculate 
on  the  order  of  magnitude  of  the  energy  exchange.  We  can  proceed  on 
this  basis  -• i  grouping  the  various  states  Into  tour  energy  levels: 
the  ground  state,  the  45  levels  or  critical  potential  levels  (some 
of  which  ere  metastable),  the  upper  excited  states,  and  the  free 
electrons.  The  processes  can  then  be  pictured  as  shown  In  Figure  4. 

Radiative  recombinations  to  all  but  the  ground  state  cause 
the  bulk  of  the  translational  and  ionisation  energy  of  the  electron 
to  be  lost,  because  the  plamaa  i»  optically  thin  for  these  wavelengths. 
Radiation  resulting  from  recombination  to  the  ground  state  will  be 
absorbed  and  reemitted  many  times  before  it  can  escape  from  the  plasma. 
This  In  effect  lengthens  the  mean  lifetime  of  the  lowest  excited  states, 
the  4S  levels,  so  that  collisions!  deexcitation  Is  the  major  means 
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of  going  from  these  I ovals  to  the  ground  state. 

Atom-atom  or  Ion-atom  suporalastlc  (deexcitation)  collisions 
In  which  tha  alactronlc  excitation  enargy  's  convartad  into 
translational  anargy  would  ba  halpful  In  reducing  tha  raqulrad  valua 
ofAF.  Howaver,  tha  procass  of  alactron-atom  collislonal  daaxcltation 
has  an  advantage  ovar  tha  haevy  partlcla  collislonal  daaxcltation, 
bacausa  of  the  higher  electron  thenal  speeds,  by  a  factor  269. 

Thus,  in  order  for  tha  two  processes  to  ba  of  equal  magnitude, 
tha  ratio  of  heavy  partlcla  density,  multiplied  by  it;;  superelastic 
cross  section,  to  the  electron  density,  multiplied  by  tha  alactron- 
atom  superelastic  cross  section,  must  ba  approximately  269. 

Some  Information  about  Inelastic  cross  sections  and 
transition  probabilities  exist.  Transition  probabilities  for  18 
transitions  between  5P-4S,  and  4P-4S  levels  ora  given  In (43), 
ranging  from  .368  x  10®  to  43.2  x  10®.  Considerable  Information 
exists  on  the  cross  section  for  production  of  excited  states  from 
tha  ground  state  (49,  50,  51).  Using  tha  principle  of  detailed 
balancing, 

(E  ♦Pt)  -  **1  E  Qct) 

1  '  / 

where  Et  *  anargy  difference  between  states  j  and  i 

•  statistical  weight  of  state*' 

in c)  •  cross  section  for  Inelastic  transfer  from  state 
4  *  to  j  caused  by  collision  of  an  electron 

of  energy  K, 

we  can  obtain  the  superelastic  collision  cross  section.  Values 

•  a  o  -|9 

range  from  about  10  '  cm.  to  10  end  lower  for  transitions 
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from  excited  states  to  the  ground  state. 

Inelastic  cross  sections  from  excited  states  to  other 

excited  states  are  difficult  to  obtain  experimentally.  None  Is 

known  to  exist  for  argon.  Atom-atom  Inelastic  cross  section  data 

are  very  scarce.  It  is  known  that  the  superelastlc  collision 

cross  section  can  be  quite  large  If  the  energy  difference  between 

the  two  states  Is  very  small,  but  as  the  difference  becomes 

large  the  cross  sections  drop  sharply(45) .  As  an  example  of 

values  for  small  energy  differences,  (46)  quotes  values  of  about 
-15  2 

10  cm.  for  transitions  In  helium  involving  energy  changes 
less  than  0, |  ev. 

These  rother  sketchy  data  lead  to  the  fol losing 
conclusions! 

(a)  Electron-atom  superelastlc  collisions  are  probably 
the  major  mechanism  of  deexcitation  from  the  critical  potential 
levels,  and  perhaps  from  the  upper  levels  to  lower  levels,  and 
to  the  ground  state. 

Superelastlc  collisions  between  heavy  particles  may 
play  a  pert  In  transitions  between  excited  states,  but  electron- 
atom  superelastlc  collisions  and  radiative  decay  are  probably 
most  Important. 

If  these  conclusions  are  valid,  most  of  the  enthalpy 
decay  of  the  electrons  Is  transferred  by  elastic  collisions  to 
the  heavy  particles.  Therefore,  we  take  Af 1  ^ 

Thus 

T e  — T *  =  v  TV 

•it 
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Values  ot  Tc-1  i  are  plotted  versus  TV  for  various  values 
of  ,n  Flflure(5),  where  the  relation  between  J-J’*  end 

has  been  given  in  Table  III. 

6.  Ionization  and  Recombination 

The  rate  of  ionization  by  electron-atom  collisions  is 
given  by  ^  H  e  ? 

jL  t  ~  J  V A  TV  Vf  ft 

where  W*  ■  number  density  of  the  neutral  atoms 

yr  •  relative  velocity  ot  electron  end  atom  'X.  \/e  - 
Tc  •  effective  lonlzution  cross  section 
we  have  assumed  the  atoms  are  stationary  coaipared  to  the  electrons. 

The  electron  distribution  can  be  approximated  by  the 
Maxwellian  distribution 


The  problem  is  therefore  that  of  simply  finding  the  correct  value 
tor  ^T(ye)  and  integrating. 

The  behavior  ot  at  high  energies  is  unimportant  to 
the  calculation  because  of  the  exponential  In  the  integral.  Petschek 
and  Byron(39)  use  the  linear  approximation 
f;  ■  7  x  |0“18  (E*  -  11.5) 
where  E'  Is  the  electron  energy  in  electron  volts. 

This  Indicates  a  threshold  at  the  critical  potential  of  argon,  and 
hence  postulates  a  multiple-step  ionization  process.  Unfortunately, 
the  data  on  which  this  expression  is  based  are  old,  and  thus 
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questionable.  More  recent  dote  ( 40,  4i)are,  untortune  ;ly,  tor 
one-step  ionization,  having  a  threshold  at  the  ionization  potential, 
but  may  be  used  to  establish-^  lower  limit  to  the  Ionization  rate: 


JjT  M  c+J'/e*'  /r.75" c  s't-lb  €,^ 

+  fti  k/0-t  £'7/6& 


Petschek  and  Byron's  cross-section  gives 


t.?r  x/c'-NiNe  jllATefJi 
,h*.  T.  -  B*./l  »  I.i*xi0r ‘K 

•  ll.5ev  a  Ij.fy/r'1 

For  equilibrium  conditions  at  I 5,000°K  and  I  atmosphere, 

(iM),  .«  x  10 22  ions/cm^-sec  (probable  value) 
For  the  higher  threshold. 


VeN* 


C'  (xl-xx,)jx 


%*  p 


4  31  m-'f}-*  K'-xtijA 

)  X ,1  /J3j(jof  •  Xu.-  ty,m*S  *  J 


re  '  n  T* 

which  results  in  a  value  tor  !5,0(XrK  and  I  atmosphere  (equilibrium) 
ot 


"  2t®  *  ,o2*  lons/cm3  -  sec  (lower  bound) 

Since  the  residence  time  of  the  gas  In  the  nozzle  used  for  the 
experiment  Is  at  least  10"*  seconds,  and  the  equilibrium  electron 
density  at  I5,00(ft<  and  I  atmosphere  would  be,  roughly,  2  x  I02*  /cm*. 
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the  existence  of  equilibrium  would  require  a  mean  ionization  rate 
of  only  2  x  I02*  ions/cm3  -sec.  Thus,  even  if  the  mean  rate  of 
Ionization  were  not  much  higher  than  that  at  the  nozzle  exhaust, 
and  we  use  the  lower  value  for  ionization  rate  f d.Ni\  It 


*  Xt'v 

appears  that  ionization  equilibrium  is  assured. 

In  the  mixing  region  we  must  consider  the  reaction  rate 
as  the  difference  between  recombination  end  lonlzatloni 

2jj£*  -  -NrV 

where  ■  "ionization"  coefficient 

0C  •  recombination  coefficient 

Using  equilibrium  values  of  Ne  and  Ns  and  taking  the 
ionization  rate  from (39),  we  obtain,  at  I5,00QPk,**  •  0.73  x  I0“*2 
cm?/ sec.  and  I  -  1 .03  x  IOM2c*.3/sec.  The  value  of  oC  so 
derived  is  in  reasonable  accord  with  the  value  expected  from  the 
assumption  of  three-body  recombination. 

Since  / JL  4  x  10*-*  /  cm  3-°K, 

~7~~~  •  A  x  I0‘3  x  ■  9  x  I02l/cm3-sec  for  the 

aLt 


extreme  cose  selected  from  (49)  . 

For  this  extreme  case,  using  (39)  and  the  above,  the  relative 
electron  densities  are  given  In  be:ow 


RELATIVE  ELECTRON  DENSITIES  FCR  EXTREME  TEMPERATURE  DECAY  RATE 
Temperature  (°K)  No/  (Ne  at  equilibrium) 


1.07 

l.ll 

1.43 

2.03 
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Tha  computed  anount  of  nonaqut I Ibrlum  alactron  danslty 
will  dacraasa  if  highar  valuas  of  tha  racomfci nation  coafficiant 
are  assunad.  In  vlaw  of  tha  uncartaintlas  inharant  in  this 
astl»ata,  it  can  only  ba  statad  that  for  tha  "axtraata"  casa  thara 
could  ba  significant  nonaqul I Ibr lum  alactron  dansity  at facts,  but 
for  tha  lass  savara  gradlants  oaasurad  in  (5)  tha  alactron 
population  is  naar  aqul I Ibriua. 

Tha  jat  is  optically  thin  in  ell  spactrel  raglons  axcapt 
thosa  of  Its  strongast  linas,  which  occur  In  tha 'far  ultravlolat. 
Hanca  sinca  tha  radiant  flux  Is  lost,  tha  nuabar  of  photo ion I  sat Ions 
will  ba  such  saallar  than  tha  nuabar  of  radiativa  racoubi nations. 

Tha  affact  of  this  labalanca  Is  to  dacraasa  tha  dagraa  of  ionisation 
(52).  Sinca  tha  thraa  body  racoafc  I  nation  Is  such  largar  than  tha 
radiativa  racoafelnatlon  this  affact  will  ba  aaall. 
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D.  Arcjet  Radiation 

Simple  total  radiation  surveys  of  the  turbulent 
argon  arcjet  were  made  using  a  collimated  radiation  probe. 
This  work  was  first  reported  in  Reference  (4).  Since  this 
time  the  results  have  been  re-examined  and  it  was  found 
that  the  total  radiation  loss  is  more  in  agreement  with 
established  theories  than  first  thought.  The  conclusion 
of  this  new  work  with  respect  to  the  turbulent  mixing 
study  is  the  same  as  previously  concluded,  i.e.,  the 
radiation  from  the  turbulent  arcjet  is  sufficiently  small 
that  it  does  not  affect  the  conclusions  of  the  turbulent 
mixing  analysis  and  experiments.  The  recent  work  did 
result,  however,  in  the  general  conclusion  that  radiation 
in  some  case  approached  ten  percent  of  the  total  turbulent 
jet  power,  substantially  higher  than  first  thought.  Also 
concluded  was  that  the  radiation  cannot  be  neglected  in 
the  mixing  studies  of  the  laminar  regime. 

No  further  discussion  of  arcjet  radiation  will  be 
made  here  since  a  report  entitled  "Analytical-Experimental 
Correlation  of  Radiation  Loss  from  an  Argon  Arcjet"  will 
be  issued  in  the  very  near  future.  The  reader  is  referred 
to  Reference  (  4  )  for  a  detailed  description  of  the  technique 
used  to  measure  radiation  power. 


IV.  CALORIMETRIC  PROBE 

A.  Introduction 

The  measurement  of  gas  properties  at  atmospheric 
pressures  and  temperatures  above  5,000°K  has  long  been  a 
serious  problem.  Devices  which  depend  on  solid-state 
properties;  e.g.,  thermocouples  or  thermometers  of  any 
type,  would  be  beyond  their  melting  points.  Simple  optical 
techniques  such  as  pyrometry  or  sodium  "D"-line  reversal 
are  limited  by  source  brightness  temperatures.  More 
advanced  spectroscopic  methods;  e.g.,  (54,  55)  etc.,  are 
not  only  delicate  and  elaborate,  but  are  subject  to  a 
number  of  diff icult-to-avoid  errors;  e.g.,  depth  of  field 
in  sources  having  high  gradients,  pressure  broadening  effects, 
difficulty  in  isolating  the  appropriate  "temperature,"  etc. 
These  errors  generally  result  in  experimental  scatter  on  the 
order  of  10  to  20%.  Other  temperature-measurement 
techniques  which  depend  on  gas  properties  such  as  electrical 
conductivity,  sound  velocity,  etc.  (56)  are  subject  to 
numerous  errors  when  applied  to  nonuniform  regions  and  have 
been  unable  to  provide  better  than  10%  repeatability. 

B.  Description 

The  probe  configuration  developed  for  the  required 
arc jet  measurements  is  diagrammed  in  Figure  6.  Construction 
of  the  probe  itself  is  of  copper,  with  stainless-steel 
supports.  Cooling  water  from  a  high-pressure  source  (up  to 
500  psi)  enters  through  the  mounting  block,  passes  up  the 
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front  stainless-steel  support,  and  through  the  outermost 
coolant  channel  to  the  probe  tip.  It  returns  via  the 
inner  coolant  channel.  Sheathed,  ungrounded  thermocouple 
junctions  are  located  precisely  at  the  probe  coolant  channel 
inlet  and  outlet. 

The  central  tube  carries  a  steady  flow  of  sample 
gas  from  the  probe  tip  past  a  thermocouple  junction  located 
precisely  opposite  the  "water  out"  thermocouple,  and  then 
through  a  gas  sample  tube  to  one  or  more  instruments  as 
described  later. 

Note  that  the  major  diameter  of  the  sampling  tube 
is  approximately  1/8  inch.  Calibration  results  on  earlier 
1/4-inch  models,  to  be  discussed  below,  clearly  indicate 
the  superior  accuracy  of  the  smaller  probe  for  local  measure¬ 
ments.  A  photograph  of  the  1/ 8-inch  probe,  also  showing  the 
exit  of  the  arc jet  used  to  provide  the  hot  gas  source,  appears 
in  Figure  7. 

C.  Methods 

1.  Temperature  Measurement 

The  calorimetric  method  used  to  determine 
gas  temperature  depends  heavily  on  a  unique  "tare"  measure¬ 
ment,  which  effectively  eliminates  errors  due  to  external 
cooling  requirements.  A  valve  in  the  gas  sample  line  is 
closed,  thus  preventing  gas  from  entering  the  probe,  and 
observations  of  coolant  temperature  rise  and  flow  rate  are 
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made.  The  valve  is  then  opened,  allowing  a  gas  sample  to 
flow  through  the  probe,  and  the  same  measurements  are 
repeated,  together  with  those  of  the  steady  gas  sample 
temperature  at  the  probe  exit  and  steady  gas  sample  flow 
rate.  The  rate  of  heat  removal  from  the  gas  sample  is 
thus  given  by  the  difference  between  the  two  coolant  a 
rates : 

*g  <hlg  “  h2g5  =  (*c  Cpc  A  Tc)flow  "  (*c  Cpc  A  Tc)no  flow 


where  w 

g 


Atc 


gas  sample  mass  flow  rate 
coolant  water  mass  flow  rate 
unknown  gas  enthalpy  at  probe  entrance 
gas  enthalpy  at  probe  exit  thermocouple 
coolant  specific  heat 

coolant  temperature  rise  =  tTc)Qut  -  (Tc>^n 


The  effectiveness  of  the  "tare"  measurement  (i.e. , 
the  heat  removal  rate  with  the  gas  sample  flow  shut  off)  is 
dependent  on  the  duplication  of  flow  conditions  near  the  probe 
tip  in  the  "flow"  and  "no-flow"  cases.  The  calibration  data 
described  below  appear  to  indicate  that  satisfactory  duplication 
is  achieved  for  the  existing  probe  geometry. 

The  unknown  gas  enthalpy  h^at  the  probe  tip  is  now 
uniquely  determined,  provided  the  gas  sample  flow  rate  v;  and 

o 

the  gas  enthalpy  hj^  at  the  probe  exit  are  known.  Since  the 
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sample  is  cooled  to  rather  low  temperatures  by  the  time  it 
reaches  the  probe  exit  thermocouple  (e.g.,  ^ 3Q0°C) ,  the 
gas  enthalpy  is  a  direct  function  of  the  temperature,  measured 
by  the  thermocouple,  and  the  gas  sample  composition,  determined 
as  discussed  in  the  next  section.  The  flow  rate  of  the  gas 
sample  may  be  measured  by  any  convenient  means;  in  the  present 
study  a  choked  orifice  was  used,  for  which  knowledge  of  the 
gas  composition  was  sufficient  to  provide  the  necessary  flow 
rate  data  (see  Figure  8). 

2.  Gas  Composition 

The  gas  sample  passes  from  the  probe  through 
a  constant-temperature  oil  bath,  and  composition  can  thus  be 
determined  at  known  equilibrium  conditions  by  various  methods. 
In  the  case  of  the  two-component  mixtures  of  interest  in  the 
present  program  (argon-helium),  a  simple  thermal  conductivity 
cell  provided  accurate  composition  data  under  steady-state^ 
gas  sample  flow  conditions  (see  Figure  8).  For  more  complex 
mixtures,  conventional  chromatography  techniques  may  be  used 
to  measure  the  necessary  composition  data. 

3.  Velocity 

At  the  time  the  "tare"  calorimeter  measurement 
is  made,  the  dead-ended  gas  sample  passage  is  filled  with 
gas  at  the  stagnation  pressure  of  the  probe-tip  location  in 
the  arc jet  stream.  A  conventional  pressure  transducer  (see 
Figure  8),  may  be  used  to  determine  this  pressure.  For  low 
subsonic  gas  flows,  the  Bernoulli  equation  provides  a  simple 
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velocity  determination: 


W-j 


l  A  p 


where  V 


lg 

/Ig 


free-stream  velocity  at  the  probe  tip 
gas  density  at  the  probe  tip 


=  measured  pressure  rise  (stagnation  pressure- 
ambient  jet  pressure) 

Note  that  in  supersonic  flow,  a  knowledge  of  the  gas 
composition  is  necessary  in  order  to  determine  the  free-stream 
Mach  number,  since  use  of  the  ideal-gas  Rayleigh  pitot  formula 
is  not  valid  for  the  partly-ionized  gas  jet.  However,  in  this 
case  the  stagnation  pressure  determination  is  essentially 
unaffected  by  ion  recombination  in  the  detached  shock  in  front 
of  the  probe,  as  was  demonstrated  experimentally  for  reacting 
flows  in  (57),  and  thus  the  probe  could  conceivably  be  useful 
in  supersonic  flows  as  well  as  for  subsonic  velocities. 

D.  Calibration 

1.  Energy  Balance 

The  fundamental  calibration  technique  is  that 
of  comparing  the  known  total  energy  of  the  arcjet  gas  with 
an  integrated  probe  energy  survey  across  the  arcjet  exit 
plane.  The  net  rate  of  energy  output  P&  of  the  arcjet  is 


P  =  E  I 
a  a  a 


&  c  (/\  T) 
a  pa  ^  ■ 
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where 

a 

=  arc jet  voltage 

*a 

=  arc jet  current 

a 

=  arc jet  coolant 

mass  flow  rate 

c 

Pa 

=  specific  heat  of  arcjet  coolant  water 

<AT>a 

=  arcjet  coolant 

temperature  rise 

account  is 

taken  of  radiation, 

since  measurements  (see 

Section  III-D)  have  indicated  that  radiation  is  usually  small 
and  reaches  about  10  per  cent  of  the  jet  energy  only  at  the 
highest  jet  powers  studied.  The  exit-plane  energy 
survey  is  made  by  using  the  probe  to  measure  gas  velocity, 
composition,  and  enthalpy  at  each  of  approximately  15 
locations  across  a  diameter  of  the  arc jet  exit  plane,  and 
determining  the  value  of  the  jet  energy  summation 


/>,  V,  h,  v  u  r 

TT  fig  lg  lg 


where  R 
r 


radius  of  the  arc jet  exit  (3/8") 
radial  coordinate 

gas  density  from  probe  measurements 
gas  velocity  from  probe  measurements 
gas  enthalpy  from  probe  measurements 
total  energy  per  unit  time 


Two  assumptions  are  needed  to  evaluate  this 
summation.  First,  a  knowledge  of  the  equation  of  state  of 
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the  partly-ionized  free-stream  gas  is  required.  It  was 
established  (2)  that  for  the  low  subsonic  velocities  of  the 
present  program,  equilibrium  of  partly-ionized  argon  at  low 
subsonic  speeds  and  one  atmosphere  pressure  was  assured.  Thus 
the  enthalpy  measurement  defines  a  temperature  which,  with  the 
measured  argon-helium  concentration,  provides  a  value  for  gas 
density  p  . 

The  second  assumption  needed  is  that  of  axial 
symmetry  of  the  arcjet,  since  a  single  diametral  survey  is 
used  to  represent  the  entire  exit  plane.  This  was  found  not 
to  be  too  stringent  a  requirement,  and  was  adequately  satisfied 
by  the  available  arcjet  source.  Note  that  some  assistance 
in  this  regard  is  obtained  by  the  use  of  a  survey  which  is 
summed  over  the  entire  diameter  instead  of  merely  a  radius, 
as  indicated  by  the  above  formulation. 

The  resultant  ratio  P^/P^,  which  represents  the 
ratio  of  integrated  probe-measured  energy  in  the  arcjet  exit 
plane  to  actual  net  energy  input  to  the  gas,  is  plotted  in 
Figure  9  for  a  number  of  arcjet  flow  rates  and  for  two  probe 
diameters.  Agreement  is  consistently  within  5%  and  the 
average  falls  within  0.5%.  Standard  deviation  for  thee 
1/8"  probe  was  calculated  to  be  3%  and  that  of  the  1/4"  probe 
was  10%.  Temperature  in  the  center  of  the  arcjet  exhaust 
exceeded  14,000°K  at  the  higher  power  levels  shown  in  this 
figure,  and  in  all  cases  the  jet  boundary  (3/8"  away)  was 
pure  helium  at  approximately  300°K. 


65 


A  second  important  calibration  reference  was  the 
mass  flow  rate  of  argon  gas  through  the  arc  generator,  as 
measured  by  a  conventional  flowmeter.  Again  using  the  probe 
exit-plane  survey  data,  the  mass-balance  ratio 

£  -it  /’Ig  Vlg  r</V*a 

-R 

was  calculated  and  is  plotted  for  a  number  of  flow  rates 
in  Figure  10.  Agreement  was  found  to  be  comparable  with 
that  of  the  energy  balance,  further  verifying  the  accuracy 
of  the  probe  technique. 

E.  Conclusions 

1.  Successful  calibration  of  a  water-cooled  gas 
sampling  probe  capable  of  steady-state  operation  at  one 
atmosphere  and  14,000°K  was  accomplished  by  energy  balance 
and  mass  balance  measurements  across  an  arcjet. 

2.  Statistical  analysis  of  a  series  of  exit-plane 
surveys  using  the  probe  resulted  in  a  standard  deviation  of 
0.03  for  a  1/8"  probe  and  0.10  for  a  1/4"  probe.  Note  that 
these  results  give  the  total  error  of  a  15-point  survey, 
and  hence  it  is  clear  that  the  accuracy  of  the  individual 
readings  making  up  each  survey  will  exceed  these  values. 

3.  Superior  scatter  behavior  of  the  smaller  probe 
indicates  its  better  suitability  for  local  measurements; 
however,  the  fact  that  the  average  of  1/4"  probe  calibration 
ratios  was  quite  close  to  unity  indicates  that  the  sampling 
technique  itself  is  satisfactory;  i.e.,  not  scale-dependent. 
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4.  The  fundamental  advantages  of  this  type  of 
probe  are  its  high  accuracy,  ease  and  convenience  of 
calibration  and  measurement,  and  ability  to  measure  several 
flow  properties  simultaneously.  It  is  able  to  perform 
these  measurements  under  not  only  the  extremes  of  pressure 
and  temperature  discussed  earlier,  but  also  in  the  presence 
of  temperature  gradients  up  to  40,000°K  per  inch. 

5.  Although  operation  of  the  copper  probe  would 
be  questionable  in  corrosive  gases,  the  application  of  gold 
plating,  which  is  easily  accomplished,  permits  such  operation. 

6.  Some  of  the  limitations  of  the  probe  technique 

are 

a.  Equilibrium  must  exist  in  the  gas  when 
temperature  is  to  be  measured,  since  the 
probe  measures  only  the  total  gas  enthalpy. 

b.  Transients  cannot  be  measured;  however,  note 
that  the  steady-state  mode  of  measurement 
has  been  observed  to  be  quite  effective  in 
averaging  the  variations  due  to  extreme 
turbulence  and  arc  fluctuations. 

c.  Spatial  resolution  is  limited  by  the  finite 
sampling  tube  diameter  (about  0.040"  for 
the  1/8"  probe). 

d.  Mass  or  energy  balance  calibrations  require 
approximately  axisymmetric  flow.  This  is 
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not  necessary  once  a  calibration  has  been 
established,  since  the  probe  measures 
local  gas  properties  directly. 

e.  In  order  to  obtain  temperature  calibrations, 
simultaneous  velocity  and  composition 
determinations  are  also  necessary.  Again, 
this  does  not  apply  once  the  calibration 
has  been  established. 

f.  The  probe  has  to  date  only  been  applied  to 
the  measurement  of  subsonic  gas  flows; 
however,  it  is  possible  that  equivalent 
property  measurements  in  supersonic  flows 
can  be  made  under  certain  conditions. 


V. 


TURBULENT  MIXING  EXPERIMENTS 


A.  Arcjet  Apparatus 

The  arcjet  used  in  these  tests  was  manufactured  by 
the  Thermal  Dynamics  Corporation.  The  basic  F-80  equipment 
(see  Figures  11  and  12)  with  peak  power  capacity  of  80 
kilowatts  was  modified  for  this  study  by  the  use  of  a  longer 
primary  nozzle  and  a  secondary  annular  helium  nozzle  (see 
Figures  13  and  14).  Power  was  obtained  from  a  440-volt 
shunt-wound  generator  driven  by  a  diesel  engine.  Five 
rectifiers  were  used  to  convert  to  direct  current. 

The  cathode  (see  Figure  12)  was  constructed  of 
ground  tungsten  and  was  water-cooled.  The  water-cooled 
copper  anode  also  served  as  the  primary  argon  nozzle.  The 
nozzles  employed  in  the  experimental  studies  were  of  0.750 
inch  diameter  straight  bore,  with  lengths  of  four  and  eight 
inches.  This  provided  L/D  ratios  of  5.33  and  10.67 
respectively  (see  Figure  15  and  16).  The  efficiencies  (net 
gas  power/electrical  input  power)  ranged  from  14%  for  the 
eight  inch  nozzle  to  about  26%  for  the  four-inch  nozzle. 

One  problem  of  special  interest  was  that  of  cathode 
concentricity  and  the  related  gas  swirl  angle.  In  order  to 
prevent  the  electrical  discharge  from  striking  the  same 
point  on  the  nozzle  at  all  times,  which  would  soon  result 
in  a  thermal  failure,  it  was  necessary  to  obtain  a  rotational 
arc  pattern.  This  was  accomplished  by  introducing  the  gas  in 
a  vortex  or  swirl.  Excessive  swirl  would  make  measurement  of 
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velocity  near  the  nozzle  boundaries  impossible  with  a 
pitot-type  device  such  as  the  probe.  Too  little  swirl 
would  result  in  asymmetric  profiles  or  thermal  failure 
of  the  nozzle.  The  optimum  swirl  for  this  study  was  found 
to  be  15°  from  axial  and  15°  from  radial.  The  swirl  plate 
used  to  produce  the  desired  flow  is  shown  in  Figure  12. 

B.  Physical  Measurements 

1.  Gas  Enthalpy:  This  was  measured  directly  by 
the  probe,  requiring  the  following  subsidiary  measurements: 

(a)  Inlet  and  exit  coolant  temperatures  from 
the  probe  (see  Figure  6).  These  were 
measured  by  standard  copper  -  constantan 
thermocouples  with  electrical  output 
observed  on  a  recording  potentiometer. 

(b)  Probe  exit  gas  temperature  (see  Figure  6). 
This  was  measured  with  a  chromel-alumel 
thermocouple  with  output  on  a  recording 
potentiometer. 

(c)  Probe  gas  flow  rate.  This  was  measured  with 
a  calibrated  critical  orifice  (see  Figure  8). 
Upstream  and  downstream  orifice  pressures 
were  measured  with  mercury  manometers, 

(NOTE:  Knowledge  of  the  probe  sample  gas 
composition  was  necessary  to  determine  the 
flow  rate) . 
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2.  Gas  Composition:  The  composition  of  the  gas 
was  determined  in  terms  of  the  helium-argon  fraction  by 
means  of  a  thermal  conductivity  cell  (see  Figure  8),  whose 
output  was  observed  on  a  recording  potentiometer.  The 
temperature  of  the  gas  within  the  cell,  cooled  by  passage 
through  an  oil  bath,  was  measured  with  a  mercury  thermometer. 

The  cell  was  calibrated  (1),  (2)  with  known,  prepared  gas 
samples  of  varying  argon-helium  fractions,  and  was  repeatably 
accurate  to  better  than  1%  (see  Reference  (5)). 

3.  Dynamic  Pressure:  The  dynamic  pressure  (Mach 
number  <  .1)  was  measured  by  means  of  a  dynamic  pressure 
transducer . (Figure  8)  and  observed  on  a  recording  potentiometer. 

4.  Probe  Radial  Position:  The  radial  position 
of  the  probe  was  measured  with  a  potentiometer  and  observed 
on  a  recording  potentiometer  (see  Figure  17). 

5.  Probe  Axial  Position;  The  axial  position  of 
the  probe  was  measured  with  a  simple  scale  amplified  by  a 
hair-line  magnifying  glass  (see  Figure  17). 

6 .  Arcjet  Voltage:  Measured  directly  with  a 
calibrated  voltmeter  (see  Figure  18). 

7.  Arcjet  Current:  Measured  directly  with  a 
calibrated  ammeter  (see  Figure  18). 

8.  Arcjet  Coolant  Temperature  Rise:  Measured  by 


copper-con stantan  thermocouples,  with  output  on  a  recording 
potentiometer  (see  Figure  19). 
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9.  Probe  Coolant  Flow  Rate :  The  water  flow  rate 
through  the  probe  (necessary  for  the  gas  enthalpy  measurement) 
was  measured  with  a  standard  calibrated  rotameter  (see  Figure  20). 

10.  Arc jet  Coolant  Flow  Rate:  The  arc jet  coolant 
flow  rate  was  also  measured  with  a  calibrated  rotameter  (see 
Figure  20). 

11.  Argon  Mass  Flow  Rate:  The  argon  volumetric  flow 
rate  was  measured  with  a  calibrated  rotameter  (see  Figure  20). 
Argon  pressure  at  the  rotameter  was  measured  with  a  simple 
bourdon  gauge,  and  argon  temperature  at  the  rotameter  was 
measured  by  a  mercury  thermometer.  These  data  were  sufficient 
to  calculate  the  mass  flow  rate. 

12.  Helium  Hass  Flow  Rate:  The  measurement  of  the 
flow  rate  of  the  secondary  gas  was  accomplished  in  the  same 
manner  as  that  of  the  primary  argon  (see  Figure  20). 

C.  Deduced  Physical  Quantities 

From  the  above  physical  measurements,  the  condition 
of  equilibrium,  and  a  known  equation(s)  of  state  for  the 
multicomponent  system  (see  Reference  (5)),  it  was  possible 
to  calculate  the  following  important  quantities: 

1.  Gas  Temperature 

2.  Gas  Composition 

3.  Flow  Velocity 

It  was  these  indirect  physical  quantities  which 
were  of  primary  importance  in  this  study. 
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D.  Nature  of  Experimental  Methods 

The  basic  purpose  of  the  experimental  program  was 
to  obtain  direct  local  measurements  over  the  entire  flow 
field  of  interest.  This  was  accomplished  in  the  following 
manner: 

After  determining  the  various  parameters  of  interest 
and  assuring  that  the  arc jet  was  yielding  reasonably  symmetric 
profiles  (this  was  essential  in  the  calibration  phase,  as  one 
diametral  pass  was  used  to  represent  the  exit  plane  conditions, 
and  also  of  general  interest  due  to  the  use  of  an  axially- 
symmetric  analysis),  the  primary  jet  and  mixing  region  were 
investigated  with  the  calorimetric  probe.  This  was  done  by 
placing  the  probe  in  a  fixed  axial  plane  and  then  moving  it 
radially  in  a  series  of  small  increments  through  the  entire 
jet.  At  each  location,  steady-state  readings  (i.e.,  the  time 
average  of  the  turbulent  fluctuations  (58),  whose  frequency 
was  observed  to  be  of  the  order  of  2000  cycles  per  second 
(see  Figure  (21))  of  all  the  physical  measurements  listed 
above  were  taken.  The  probe  then  (typically)  was  moved 
to  about  15  such  locations  in  the  given  axial  plane.  A 
different  axial  position  was  then  selected  and  the  procedure 
repeated.  The  number  of  such  axial  positions  was  typically 
six  or  seven. 

The  procedure  outlined  above  resulted  in  data 
for  one  fixed  set  of  flow  parameters.  The  need  for  variation 
of  such  basic  flow  parameters  as  peak  temperature, 
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peak  velocity,  velocity  ratio,  scale  and  intensity  of 
turbulence,  etc.  accounted  for  the  large  amount  of  experimental 
information  collected.  For  reasons  of  brevity  (see  Table  2) 
only  a  very  small  portion  of  the  data  obtained  has  been 
included  in  this  report. 

From  the  previously  stated  physical  measurements 
it  was  possible  to  determine,  for  various  values  of  the 
flow  parameters, 

(a)  Velocity,  temperature  and  helium 
concentration  profiles  at  various  axial 
positions  (see  Figures  22,  24  and  26). 

(b)  The  axial  decay  of  velocity,  temperature 
and  helium  concentration  (see  Figure  28,  29, 
30). 

(c)  The  axial  and  radial  gradients  of  these 
properties  over  the  flow  field  (see 
Figures  31,  32,  33). 

(d)  The  radial  spreading  of  these  same  quantities 
as  functions  of  axial  position  (see  Figure  34). 

The  profiles,  axial  centerline  decay  and  axial 
centerline  gradients  presented  no  particular  difficulties. 

However,  as  the  determination  of  the  radial  spreading  of 

the  field  variables  was  not  quite  so  simple,  a  brief  description 

of  this  matter  is  in  order. 
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Essentially,  the  determination  of  the  radial 
spreading  of  a  field  variable,  say  temperature,  stemmed  from 
inspection  of  the  radial  temperature  profile.  Since  the 
energy  or  temperature  outer  boundary  is  defined  as  the 
locus  of  points  nearest  the  centerline  where  T  =  To*  ,  then 
by  determining  the  value  of  the  radial  position  at  which 
this  condition  is  met,  for  a  profile  at  a  given  axial  position, 
the  outer  boundary  can  be  related  to  the  axial  position.  This 
was  done  for  all  the  various  axial  positions  and,  of  course, 
the  other  field  variables. 

Nonetheless,  a  major  difficulty  still  remained, 
for  the  exact  point  at  which  a  variable  "reaches"  its 
asymptotic  limit  is  very  hard  to  ascertain.  Thus  a  value 
slightly  removed  from  this  limit  was  selected,  and  by 
repeating  this  procedure  for  different  arbitrary  values  it 
was  possible  to  extrapolate  to  the  true  boundary  criteria. 

This  procedure  was  employed  for  determination  of  all 
experimental  boundary  spreading. 

E.  Results 

1.  The  integral  mixing  analysis  gives  good  agreement 
with  experimental  results  based  on; 

a.  Prediction  of  radial  boundaries  of 

concentration,  temperature,  and  velocity 
(see  Figure  34). 
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b.  Axial  decay  of  concentration,  velocity  and 
temperature  (see  Figures  28,  29,  30). 

c.  Axial  gradients  of  concentration,  velocity 
and  temperature  (see  Figures  29,32,  33). 

d.  Radial  profiles  of  concentration,  velocity 
and  temperature  (see  Figures  22,  23,  24, 

25,  26  and  27). 

2.  The  "cusped  parabaloid"  shapes  of  the  spreading 
boundaries  (Figure  34)  resulting  from  the  numerical  solution 

of  the  integral  equations  differ  from  the  conical  boundaries 
calculated  by  Squire  and  Trouncer  (27)  and  obtained  by  Warren 
(15)  at  low  temperatures,  but  are  in  qualitative  agreement  with 
the  statistical  theory  due  to  Taylor  (10). 

3.  The  Prandtl  number  and  Schmidt  number  are  indeed 
less  than  unity,  while  the  Lewis  number  is  greater  than  unity 
(see  Table  3  and  boundary  relationship  of  Figure  1). 

4.  Both  stable  and  "unstable"  analytical  solutions 
were  obtained.  "Unstable"  solutions  occur  when: 

a.  The  various  boundaries  for  spreading  are  too 
widely  separated. 

b.  The  driving  parameters  for  mass  diffusion 
momentum  transfer  or  heat  transfer  are  too 
large. 

c.  Any  one  such  driving  term  is  larger  by  an 


76 


order  of  magnitude  than  the  other  two. 

5.  The  axial  derivatives  of  all  the  field  variables 
are  linearly  proportional  to  the  axial  derivatives  of  their 
respective  boundaries  (see  Reference  (5)).  This  is  a  direct 
consequence  of  the  similarity  assumption  regarding  radial 
profiles. 

6.  The  Reichardt  hypothesis,  used  to  formulate  the 
driving  term  in  the  integral  momentum  equation,  provides  good 
agreement  with  experiment  for  this  turbulent,  highly  non- 
isoenergetic  case. 

7.  Analytically,  the  magnitude  of  radial  spreading 
was  found  to  decrease  slightly  with  increased  temperature. 

This  effect  was  far  more  pronounced  experimentally  (see  Result 
No.  12).  One  possible  explanation  for  this  behavior  could 

be  the  analytical  assumption  of  mean  values  for  the  diffusion 
coefficient  and  the  thermal  conductivity. 

8.  Concentration  was  found  to  spread  more  rapidly 
than  both  temperature  and  velocity,  and  temperature  to  spread 
more  rapidly  than  velocity,  in  excellent  agreement  with  the 
analytical  result  (see  Figure  34). 

9.  Experimental  profiles  of  concentration,  velocity 
and  temperature  in  the  mixing  region  were  observed  to  be 
strikingly  similar  (see  Figures  22,  24,  26),  supporting  the 
similarity  assumption  used  in  tie  analysis. 

10.  A  potential  core  for  concentration,  velocity,  and 
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temperature  was  observed.  Its  length  was  approximately  one 
jet  diameter  for  all  the  field  variable!,  and  was  relatively 
insensitive  to  temperature  level.  Profile!  within  this 
core  were,  of  course,  definitely  dissimilar  (see  Figures  22 
through  33). 

11.  Experimental  spreading  boundaries  showed  a  "cusped 
paraboloid"  shape  at  typical  operating  plasma  temperatures, 

and  a  conical  shape  in  the  isoenergetic  (unheated)  cases  (see 
Figure  34 ) . 

12.  The  magnitude  of  the  boundary  spreading  of  all 
three  quantities  was  observed  to  decrease  with  the  onset  of 
ionization  and  then  became  constant  at  ionization  levels  greater 
than  10%  (see  Figures  35,  36,  37,  38,  and  40,  and  Result  No.  7). 

13.  The  axial  decay  of  concentration,  velocity 

and  temperature,  while  in  excellent  qualitative  agreement  with 
other  experimental  results  (15),  (17),  (19),  (59)  and 
theoretical  calculations  (15),  (19),  and  (27),  is  much  more 
rapid  for  this  high  temperature  case  than  in  the  earlier 
isoenergetic  or  near-isoenergetic  cases  reported  in  the 
literature.  Typical  "length  to  approximate  full  decay"  in 
this  study  was  eight  radii  (see  Figures  28  ,  29,  30).  The 
above  papers  obtained  typical  lengths  of  40  radii. 

14.  The  plasma  was  highly  turbulent,  with  no  notice¬ 
able  transverse  oscillations  and  with  a  longitudinal  frequency 
of  approximately  2000  cycles/second  for  the  eight-inch  nozzle. 
There  was,  however,  a  significant  effect  of  nozzle  L/D  ratio, 
in  that  turbulence  in  the  four-inch  nozzle  was  not  only  of 
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much  greater  intensity  than  in  the  eight-inch  nozzle,  but 
also  exhibited  definite  transverse  displacements  (see  Figure 
21) . 

15.  The  three  transfer  mechanisms,  mass  diffusion, 
momentum  transfer,  and  energy  transfer  were  found  to  be 
extremely  effective  in  the  case  of  hot,  partially  ionized 
argon  mixing  in  a  turbulent  fashion  with  a  cool,  inert, 
helium  stream.  This  is  indicated  by  the  short  potential 
core  and  the  very  rapid  axial  decay  rateB  discussed  in 
Result  No.  13. 

16.  The  concentration  and  temperature  boundaries 
were  independent  of  variations  of  A ,  the  secondary  stream 
to  primary  stream  velocity  ratio  (see  Figure  39).  The 
momentum  boundary  R*  showed  a  weak  linear  variation  with  A 
in  contrast  to  the  results  of  Squire  and  Trouncer  (27).  How 
ever,  the  range  of  ?\  was  limited  and  this  trend  was  not 
conclusively  established. 
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F.  List  of  Symbols 


a 

aij 

A» 

Az 

C 

9 

h 

X 

X 

L 

M 

Wr« 

Npr 

Mse 

Nte 

^Pe 

Klia. 

P 

r 

R 

T 

U 


Radius  of  primary  nozzle 

Matrix  integral  coefficient  (21  in  all) 

Spreading  relation  parameter  (Ref.  5,  Equation  B-4b) 
Spreading  relation  parameter  (Ref.  5,  Equation  B-4b) 
Concentration  (particle  fraction  of  a  species) 

Sonic  velocity 


Function 

of 

CHe* 

8 

( see 

Ref. 

5) 

Function 

of 

CHe* 

8 

( see 

Ref. 

5) 

Function 

of 

CHe* 

8 

(see 

Ref. 

5) 

Thermal  conductivity 
Prandtl  mixing  length 
Length  of  potential  core 
Mach  number 
Reynolds  number 
Prandtl  number 
Schmidt  number 
Lewis  number 
Peclet  number 
"Jacobs"  number 

Dimensionless  number  (see  Ref.  5) 
Physical  radial  coordinate 
Dimensionless  radial  coordinate  =  r/a 
Absolute  temperature 
Mean  axial  velocity  component 
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u 

X 


o<(diyj)  = 
= 

<3  - 

r  = 

A  = 

d 
£ 

t  - 

G 

K 

* 

TT 

e  = 

X 

b  = 

V 


Dimensionless  axial  velocity  component  (see  Ref.  5) 
Physical  axial  coordinate 

Greek  Letters 


Function  of  matrix  coefficients  defined  by 
Equation  B-55,  Ref.  5 

Function  of  matrix  coefficients  defined  by 
Equation  B-55,  Ref.  5 

Dimensionless  parameter  defined  in  Ref.  5 

Dimensionless  parameter  defined  in  Ref.  5 

Denotes  finite  difference 

Denotes  partial  differentiation 

Turbulent  momentum  exchange  coefficient 

Dimensionless  axial  coordinate  =  x/a 

Dimensionless  temperature  (see  Ref.  5) 

Dimensionless  turbulent  exchange  constant 
(see  Ref.  5) 

Dimensionless  velocity  ratio  (see  Ref.  5) 

3.14159 . 

Density 

Turbulent  shear  stress 
Laminar  viscosity 
Kinematic  viscosity 


Subscripts 


Outer  boundary 
Helium 

Turbulent  value  of 
For  diffusion 

Denotes  typical  ith  row,  jth  column  matrix  element 
Values  in  free  outer  stream 

Superscripts 

Concentration 

Momentum 


Energy 
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TABLE  I 

kl£T.£f.  EXPERIMENTAL  Obit  RUNS 


CASE 

TORCH 

CURRENT 

AMPERES 

ARGON 

FLOW 

S.C.F.H. 

HELIUM 

FLOW 

S.C.F.H. 

NOZZLE 

LENGTH 

INCHES 

PEAK 

TEMPERATURE 

°R 

PEAK' 

VELOCITY 

♦t/sec 

PEAK 

ENTHALPY 

bfu/lb 

PEAK  * 
PERCENT 
IONIZATION 

D-4 

1500 

350 

5.0 

8.00 

22,200 

440 

5520 

i3.3 

D-5 

1750 

350 

5.0 

8.00 

23,000 

465 

6600 

25.0 

0-6 

1250 

350 

5.0 

8.00 

.22,000 

400 

-5300 

16.2 

D-7 

1000 

350 

5.0 

8.00 

21,900 

370 

5200 

15.8 

D-B 

'  5.0 

350 

5.0 

8.00 

20,300 

280 

3700 

7.2 

D-17 

0 

350 

5.0 

8.00 

530 

30 

0 

0 

D-ll 

1500 

350 

10.0 

8.00 

22,200 

390 

5530 

13.3 

D-12 

1500 

350 

15.0 

8.00 

21,500 

400 

4750 

10.0 

D-14 

1500 

350 

20.0 

8.00 

21,400 

440 

4650 

12.5 

0-10 

1500 

450 

5.0 

8.00 

22,800 

480 

6000 

22.0 

C-9 

1500 

400 

5.0 

8.00 

23,300 

520 

7000 

28.0 

D-15 

1500 

320 

5.0 

8.00 

23,800 

600 

7800 

34.6 

D-16 

0 

350 

5.0 

4.00 

530 

30 

0 

0 

0-18 

500 

350 

5.0 

4.00 

21,800 

360 

5010 

11.5 

0-3 

1000 

350 

5.0 

4.00 

17,000 

680 

2300 

00.5 

0-20 

250 

350 

5.0 

8.00 

20,000 

230 

3500 

5.3 

0-21 

2000 

■  350 

5.0 

8.00 

*  Note:  Refers  to  Cjo ,o 
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TABLE  2 

usuEmnm  cases  comsioebeo 


CASE 

TOjO  (°R) 

NR«D 

Ja 

*1 

a2 

X 

- 

\»,o 

L 

NN 

(MESH  SIZE) 

1 

18,000 

500 

2000 

.0005 

0.0120 

0.0150 

nra 

n 

0.01 

2 

18,000 

500 

1000 

.0005 

0.0120 

0.0150 

EE 

mm 

0.01 

3 

18,000 

500 

1000 

.0010 

0.0120 

0.0150 

0.05 

0.05 

S9 

0.01 

4 

18,000 

1000 

2000 

.0005 

0.0120 

0.0150 

0.05 

0.05 

mm 

0.01 

5 

18,000 

2000 

4000 

.0002 

0.0120 

0.0150 

0.05 

0.05 

0 

o.ot 

6 

18,000 

500 

1000 

.0001 

0.0000 

0.0000 

0.05 

EEH 

0 

0.01 

7 

18,000 

500 

1000 

.0001 

0.0040 

0.0060 

0.05 

EE1 

0 

0.01 

8 

18,000 

500 

iOOO 

.0010 

0.0060 

0.0080 

0.05 

0.05 

0.00 

0.01 

9 

18,000 

500 

1000 

.0010 

0.0070 

0.0100 

0.05 

0.00 

0.01 

10 

18,000 

500 

1000 

.0010 

0.0060 

0.05 

0.05 

0.00 

0.04 

II 

18,000. 

500 

1000 

.0010 

0.0040 

0.0060 

0.10 

0.03 

0.20 

0.10 

12 

18,000 

500 

1000 

.0010 

0.0040 

0.0060 

0. 15 

0.03 

0,50 

0.10 

13 

18,000 

500 

1000 

.0010 

0.0040 

0.0060 

0.20 

0.05 

1.00 

0.10 

14 

18,000 

500 

1000 

.0010 

0.0040 

0.0060 

0.05 

0.05 

2.00 

0.01 

15 

15,000 

500 

500 

.0025 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

16 

20,000 

500 

500 

.0025 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

17 

25,000 

500 

500 

.0025 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

18 

21,600 

400 

400 

.0030 

0.0050 

0.0070 

0„05 

0.05 

2.00 

0.01 

19 

21,600 

300 

300 

.0040 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

20 

.  _ 

21,600 

250 

250 

.0040 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

21 

21,600 

200 

200 

.0050 

0.0050 

0.0070 

BEEi 

2.00 

0.01 

22 

21,600 

350 

300 

.0060 

0.0050 

0.0070 

EG  3 

EK3 

2.00 

0.01 

23 

21,600 

350 

300 

.0050 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

24 

21,600 

350 

300 

.0040 

0.0050 

0.0070 

0.05 

0.05 

2.00 

0.01 

25 

14,300 

350 

300 

.0040 

0.0050 

0.05 

0.05 

Z.00 

0.01 

26 

9,000 

350 

300 

.0040 

0.0050 

0.0070 

0.05 

0.06 

2.00 

0.01 

27 

7,200 

350 

300 

.0040 

0.0050 

0.0070 

0.05 

0.07 

2.00 

0.01 

28 

3,600 

350 

300 

.0040 

0.0050 

0.0070 

0.05 

0.14 

2.00 

0.01 

29 

1,800 

350 

300 

.0040 

0.0050 

0.0070 

0.05 

0.28 

2.00 

0.01 

TABLE  3 


VALUES  OF  DIMENSIONLESS  PARAMETERS  FOR  A  TYPICAL  ANALYTICAL  CASE* 


fMYSICAL  MEANING 

QfflhmflN 

sasoL 

VALUE 

Ray no Ids  Number 

■ 

Inertial 

Z  V.,.  <»■  P 

9 

ffee 

•  BIO 

Viscous 

h 

Prendtl  Number 

m 

Haaanfii  Trcntfr 

-A-h- 

m 

Npr 

•  0.82 

Thermal  Transfer 

jk 

Schmidt  Number 

m 

Momentun  Transfer 

—A — 

9 

NSc 

•  0.74 

. 

Mass  Transfer 

PB 

Lewis  Number 

m 

Hfl»a  Tranafr 

+i£ 

9 

NLe 

-  1.11 

Thermal  Transfer 

' 

Pec  let  N unbar 

m 

^"Pr 

9 

NPe 

•  867 

"Jacobs"  Number 

9 

^.“sc 

9 

NJn 

-  600 

Turbulent  Reynolds  • 

2  u.  r„" 

Number 

Viscous 

€ 

K. 

-  TOO 

<» 

m 

P 

-  0.0005 

Thermal  Energy 

o 

V 

fro 

e 

m 

-1 

Mach  Number 

9 

-I*iflfcLt* _ 

Sonic  Velocity 

Co^e 

• 

-  0.10 

Analytical  Casa  24 
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TABLE  4 

jttfc  CHARACTERISTIC  EXPERIMENTAL  VALUES 
Gas  Properties: 


1) 

Psek  Temperatures 

14,000  °K  »  25,200  °R 

2) 

Peak  Enthalphyj 

5,800  cel/gram  *10,500  btu/lb 

3) 

Density: 

40  gram/m3  »  2  x  I0"3  lb/ft3 

4) 

%  Ionization  (Paek)i 

34* 

5) 

Paak  Primary  Valocifyt 

215  m/ sac  *  700  ft/ sac 

6) 

Paak  Mach  Number: 

0.10 

7) 

Paak  Secondary  Velocity: 

22  m/sac  •  72  ft/sac 

8)  Primary  Argon  Flow  Rata: 

Torch  Propart  s' ssi 

2  grams/sac  •  4  a  IQ"3  lb/ sac 

9) 

Paak  Currant) 

2000  amperes 

10) 

Peak  Voltagai 

40  volts 

II) 

Paak  Nat  Powers 

20  K.W.  <25*  efficiency) 

12) 

Torch  Coolant  Flows 

1  kg/sac  •  2.2  Ib/sec 

13)  Torch  water  Temperature  Rises 

Proba  Properties: 

30  °C  -  54  °R 

14) 

Paak  Probe  Exit  Gas  Tamparetura: 

600  °K  -  1080  °R  «r  650  °F 

15) 

Proba  Sample  Gas  Flow  Rotes 

3  x  I0~2  9/ sac  *7*  I0"5  |b/sec 

16) 

Coolant  Tamparetura  Risat 

30  °K  ■  54  °R 

17) 

Average  Gas  Velocity  in  Sample 

100  a/ sac  ~  30O  ft/sac 

18) 

Coolant  Flow  Velocity) 

7  m/sec  «—  20  ft/sac 

19) 

Proba  Coolant  Flow  Rates 

15  gram/ sac  *  J*  IQ"2  lb/ sec 

20) 

Average  Proba  Mall  Heat  Transfer  Rates 

1  btu/in2-see 

21) 

Typical  Copper  Wall  Tamparature) 

360  «  650  °R 
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CALORIMETRIC  PROBE  INSTALLATION 
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DIAGRAM  OF  F-80  ARCJET  TORCH 
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FIGURE  12 


F-80  ARCJET  TORCH,  SWIRL  PLATE,  AND  CATHODE 
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DIAGRAM  OF  4"  AND  8"  ARGON  NOZZLES 

FIGURE  15 
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AND  8  ARGON  NOZZLES 


FIGURE  17 


PROBE  DRIVE  MECHANISM 
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INSTRUMENT  PANEL 
AND  EXPERIMENTAL  APPARATUS 

FIGURE  20 
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ALYTICAL) 


DIMENSIONLESS  TEMPERATURE 
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FIGURE  26 


TYPICAL  RADIAL  TEMPERATURE  PROFILES  (ANALYTICAL) 
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TYPICAL  AXIAL  DECAY  OF  VELOCITY 
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TYPICAL  AXIAL  DECAY  OF  CONCENTRATION 
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FIGURE  31 


TYPICAL  AXIAL  GRADIENTS  OF  VELOCITY 
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FIGURE  36 


TYPICAL  RADIAL  SPREADING  OF  VELOCITY,  CONCENTRATION 
TEMPERATURE  WITH  RESPECT  TO  PEAK  ENTHALPY 
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Mass.  Inst,  of  Tech. 

Department  of  Mechanical  F.ng'g 
Cambridge  39,  Mass. 

Attn:  Engineering  Library  (I) 

University  of  Michigan 
228  W.  Engineering  Bui Iding 
Ann  Arbor,  Michigan 

Attn:  Eng'g  Library  (I) 

University  of  Minnesota 
Mechanical  Engineering  Dept. 
Minneapolis,  Minnesota 
Attn:  Library  (I) 

Union  College 

Department  of  Mechanical  Eng'g 
Schenectady  8,  New  York 
Attn:  Dr.  G.  M.  Ketchum 
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University  of  New  Mexico 
Department  of  Mech.  Engineering 
A I buqueroue,  New  Mexico 
Attn:  Library  (I) 

Prof.  Richard  E.  Balzh  $er 
University  of  Michigan 
Dept,  of  Chem  4  Metallurgical  Eng'g 
Ann  Arbor,  Michigan  (I) 

New  York  University 
College  of  Engineering 
Dept,  of  Mechanical  Engineering 
University  Heights 
New  York  53,  New  York 
Attn:  Eng'g  Library  (I) 

> 

Northwestern  Uni  vers I  tv 
Mechanical  Engineering  Department 
Evanston,  Illinois 
Attn:  Library  (I) 

Project  SQUID 
University  of  Virginia 
Charlottesville,  Virginia 

Attn:  Dr.  J.  Scott,  Dir.  (I) 

Purdue  University 
School  of  Aero.  4  Eng'g  Sciences 
Lafayette  ,  Indiana 
Attn;  Library  (I) 

University  of  Santa  Clara 
Santa  Clara,  California 
Attn:  Library  (I) 

Stanford  University 
Stanford,  California 
Attn;  Eng'g  Library  (I) 

Syracuse  University 
Syracuse,  no*  York 
Attn;  Eng'g  Library  (I) 

Tufts  University 

Mechanical  Engineering  Depa  tment 
Medford  55,  Mass. 

A+tn:  Eng'g  Library 

Un' versify  of  Wa&.,.ngtoi 
College  of  Engineering 
Seattle  5,  Washington 
Attn:  Eng'g  Library  (■; 
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University  of  Wisconsin 
Department  of  Mechanical  Eng'g 
Madison,  Wisconsin 

Attn:  Eng'g  Library  (I) 

Aerochem  Research  Laboratories 
Inc. 

Princeton,  New  Jersey 
Attn:  Dr.  H.  F.  Calcote,  Pres.  (I) 

Aerojet-General  Corp. 

Space  Technology  Division 
Azusa,  California 
Atrn:  Library  (I) 

Aero Jet -Genera I  Corp. 

Liquid  Rocket  Plant 
Sacramento,  California 
Attn:  J.  Moise  (I) 

Aerojet-General  Nucleonics 
Fostotla  Way 
San  Ramon,  California 
Attn:  Library  (|) 

Aeronutronic  Systems,  Inc. 

Newport  Beach,  California 

Attn;  Dr.  T.  A.  Bergstrol I  (I) 

Aerophysics  Development  Corp. 

P.  0.  Box  567 

Pacific  Palisades,  California 
Attn:  Library  (I) 

Air  Reduction  Companv 
Murray  Hill,  New  Jersey 
At+n:  Library  (I) 

A I Research  Mfg.  Company 
9851-9951  Sepulveda  B I vd. 

Los  Angeles  45,  California 
Attn:  Library  (l) 

A I  I is-Chaimers  Mfg.  Company 
General  Machinery  Division 
Milwaukee  I,  Wisconsin 
Attn:  Mr.  R.  C.  Allen 

Dir.  of  Mech.  Eng'g  ( I) 

Argonne  National  Laboratories 
P.  0.  Box  299 
Lemont,  Illinois 
Attn:  Technical  Library  (I) 
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Amer Ican-Standard 
Atomic  Energy  Division 
369  Whisman  Road 
Sunnyvale,  California 
Attn:  Hr.  Neal  F.  Lansing  (I) 

Arnold  Equipment  Corporation 
3080  Main  Street 
Buffalo  14,  New  York 

Attn:  Mr.  R.  G.  Tessmer  (I) 

Aro,  Inc. 

Technical  Information  Branch 
Tullahoma,  Tennessee 

Attn;  Mr.  G.  E.  Randall 

Atomic  Power  Dev.  Ass.  Inc. 

191 1  First  Street 
Detroit  26,  Michigan 
At+n;  Library  (I) 

AVCO  Research  Laboratory 
Wi Imington,  Massachusetts 
Attn;  Dr.  R.  R.  John  (I) 

AVCO -Everett  Research  Laboratory 
2385  Revere  Beach  Parkway 
Everett  49,  Massachusetts 
Attn:  Technical  Library  (I) 

Dr.  Arnold  Goldburg  (I) 

The  Babcock  and  Wilcos  Company 
Research  and  Development  Center 
P.  0.  Box  835 
A 1 1  lance,  Ohio  ( I ) 

Bell  Aircraft  Corporation 
Space  Flight  and  Missiles  Blv. 
Buffalo,  New  York 
Attn;  Library  (I) 
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California  Research  Corporation 
Richmond,  California 

Attn:  Dr.  J.  H.  Macpherson  <l> 

Carrier  Corporation 
300  South  Geddes  Street 
Syracuse  I,  New  York 
Attn:  Engineering  Library 


Cornell  Aer nautical  Lab. 

4455  Genesee  Street 
Buffalo  21,  New  York  (I) 

Chrysler  Corporation 
Automotive  Research,  Dept.  981 
12843  Greenfield  Avenue 
Detroit  27,  Michigan 
Attn;  Mr.  A.  Charlamb 
Research  Eng.  (I) 

Detroit  Edison  Company 
Nuclear  Power  Division 
Detroit,  Michigan  (I) 

Eitel-McCul lough,  Inc. 

798  San  Mateo  Avenue 
San  Bruno,  California 
Attn:  Mr.  J.  S.  McCullough 
Mgr.  Res.  (I) 

Electric  Boat  Company 
Groton,  Conn. 

Attn:  Mr.  E.  S.  Dennison  (I) 

Electro-Optical  Systems, Inc. 

3016  E.  Foothill  Bivd. 

Pasadena,  California 
Attn;  Library  (I) 

Fairchild  Engine  and  Airplane 
Corporation 

Fairchild  Aston Ics  Division 
Wyandanch,  Long  Island,  N.  Y. 

Attn:  Jbrary  (I) 

Foster  wheeler  Corporation 
165  Broadway 
New  York  6,  New  York 
Attn:  Mr.  E.  L.  Daman 

Dir.  Res.  and  Eng'g  (I) 

General  Electric  Company 
Atomic  Energy  Division 
San  Jose,  California 
Attn;  Mr.  N.  E.  Tippets  (I) 

General  E'ectrlc  Company 
Cincinnati  15,  Ohio 
Attn:  Aircraft  Gas  Turbine  Dlv. 

Central  Horary  Build.  305 
(I) 
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Ramo  Wooldrige  Division 
8433  Fa  I  I  brook  Avenue 
Canoga  Park,  California 

Attn:  Dr,  David  B.  Langmuir  (I) 

Stanford  Research  Institute 
Menlo  Park,  California 

Attn:  Technical  Library  (I) 

Sundstrand  Denver 
2480  W.  70th  Avenue 
Denver  21,  Colorado  (I) 

Thermal  Dynamics  Corporation 
Hanover,  New  Hampshire 

Attn:  Mr.  J.  Browning,  Pres.  (I) 

Thiokol  Chemical  Corporation 
Elkton  Division 
Elkton,  Maryland 
Attn:  Dr.  Donald  D.  Thomas  (I) 

Thompson  Ramo-Wooldr Idge,  Inc. 
Research  and  Development 
TAPCO  Division 
Wi I loughDy,  Ohio 

Attn:  Mr.  W.  A.  Compton,  Mgr.  (I) 
Mr.  V.  P.  Kovacik  (I) 

Umon  Carbide  Nuclear  Company 
Oak  Ridge  Gaseous  Diffusion  Plant 
Plant  Records  Department 
P.  0.  Box  P 

Oak  Ridge,  Tennessee  (I) 

United  Aircraft  Corporation 
The  Library 
400  Ma i n  Street 
East  Hartford  8,  Conn.  (I)  • 

Vitro  Laboratories 
200  Pleasant  Valley  Way 
West  Orange,  New  Jersey 

Attn:  Dr.  C.  Dudley  Fitz  (I) 

Westinghouse  Electric  Corporation 
Atomic  Power  Division 
P.  0.  Box  1468 
Pittsburgh  30,  Pa. 

Attn:  Tech.  Library  (I) 

Westinghouse  Electric  Corp 
Research  Laboratory 
East  Pittsburgh,  Pa. 

Attn:  Dr.  Stewart  Way  (I) 


Stewart-Warner  Corporation 
1514  Drover  Street 
Indianapolis  7,  Indiana 
Attn:  Mr.  R.  D.  Randall 

Mgr.  Research  (I) 

Nucor  Research,  Inc. 

2421  Wolcott  Avenue 
Ferndale  20,  Michigan 

Attn:  Mr.  A.  R.  Bobrowekl 
V.  Pres.  (I) 

Heat  Division 
Mech.  Eng.  Res.  Lab. 

D.S.I.R.,  East  Kilbride 
near  Glasgow,  Scotland 
Attn:  Mr.  A.  J.  Ede 
SEND  TO 

Office  of  the  Asst.  Neva I 
Attache  for  Research 
Naval  Attache,  American  Embassy 
Navy  No.  100,  Fleet  P.  0. 

New  York,  New  York  (I) 

Dept,  of  Mechanical  Engineering 
Queen  Mary  Col  lege 
London  E.  I.,  England 
Attn:  Mr.  E.  J.  LeFevre  (I) 

Minneapol is-Honeywel I  Regulator 
Company 

Research  Center 

500  Washington  Avenue  South 

Hopking,  Minnesota 

Attn:  Mr.  A.  J.  Oberg  (I) 

Ingersol 1-Rand  Company 
Painted  Post,  New  York 
Attn:  John  B.  Caret, 
Librarian  (I) 

Dr.  Clerk  A.  Dunn 
Director,  Eng'g  Research 
Oklahoma  State  University 
Still  water ,  0k  I ahoma  ( I ) 

ASTRO 

A  Div.  of  the  Marquardt  Corp, 
16555  Sat  I coy  Street 
Van  Nuys,  Cal  ifornia 
Attn:  Edward  T.  Pitkin,  Mgr. 
Space  Propulsion  Dept. 

Dr.  Julio  Cordero,  AVCO  RAD 
Wilmington  Massachusetts  (I) 
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General  Electric  Company 
General  Engineering  Laboratory 
Schenectady  5,  New  York 
Attn:  Mrs.  Florence  F.  Buckland  (I) 
Or.  H.  T.  Nagamatsu,  Dir., 
Hypersonic  Lab.  ( I ) 

General  Metals  Corporation 
Enterprise  Division 
18th  and  Florida  Streets 
San  Francisco,  California 

General  Motors  Corporation 
Alii  son  Division 
Indianapolis  6,  Indiana 
Attn:  Mr.  R.  M.  Hazen 
Dir.  of  Eng.  (I) 

General  Motors  Corporation 
Research  Labs,  01  vis  ion 
Detroit  2,  Michigan 
Attn:  Library  (I) 

Hughes  Aircraft  Company 
Florence  and  Teake  Streets 
Culver  City,  California 
Attn:  Tech.  Library  (I) 

Lawrence  Radiation  Laboratory 
Box  808 

Livermore,  California 
Attn:  Dr.  T.  Merkle  (I) 

Linde  Company 

Div.  of  Union  Carbide  Corp. 

Speedway  Laboratories 
1500  Pol co  Street 
Indianapolis  25,  Indiana 
Attn;  Librarian  (I) 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division 
Aerodynamics  Department 
Sunnyvale,  California 
Attn:  Mr.  G.  A.  Etema  -  (I) 

Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attn;  Mr.  W.  Corliss, 

Nuclear  Division  ( I ) 


MHD  Research 

1571  Placentia  Avene. 

Newport  Beach,  California 
Attn:  Or.  V.  Blackman  (I) 

North  American  Aviation,  Inc. 
International  Airport 
Los  Angeles  45,  California 
Attn:  Mr.  M.  S.  Sulkln  ( I) 

Plasmadyne  Corporation 
3639  South  Main  Street 
Santa  Ana,  California 
Attn:  Library  (I) 

Radio-Corporation  of  America 
Astro-Electronic  Products  Div. 
Route  I 

Princeton,  New  Jersey 
Attn:  Mr.  $.  Sternberg 
Chief  Engineer  (I) 

Reaction  Motors  Division 
Thlokol  Chemical  Corporation 
Denville,  New  Jersey 
Attn:  Dr.  E.  Seymour 

Dir.  Research  (I) 

Plasma  Propu.sion  Laboratory 
Republic  Avaitlon  Corporation 
Farmingdale,  Long  Island,  N.  Y. 
Attn:  Mr.  A.  E.  Kunen,  Mgr.  (I) 

Rocket dyne  Division 
6633  Canoga  Avenue 
Canoga  Park,  Cali  for nJa 
Attn:  Librarian  ( I ) 

Dr.  R.  Dl I laway, 

Nuc I  ionics  Div.  ( I ) 

Shell  Development  Company 
4560  Horten  Street 
Emeryville,  California 
Attn:  Dr.  C.  R.  Oarbett  (I) 

Solar  Aircraft  Company 
San  Diego  12,  California 
Attn;  Mr.  P,  A.  Pitt, 

Chief  Engineer  (I) 
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U.  S.  Atomic  Energy  Commission 
Tech.  Into  Service  Extension 
P.  0.  Box  62 
Oak  Ridge,  Tennessee  (1) 

■Y-iJi 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
Navy  No.  100,  Fleet  P.  0. 

New  York,  New  York 

Arthur  D.  Little,  Inc. 

15  Acorn  Park 
Cambridge  40,  Mass. 

Attn:  Miss  Dorothy  E.  Hart 
Librarian  (I) 

Conva i r 

A  Di v.  of  General  Dynamics  Corp. 

Fort  v/orth  I,  Texas 

Attn:  Div.  Research  Library  (I) 

Aerojet-General  Corporation 
P.  0.  Box  296 
Azusa,  California 
Attn:  Robert  D.  Waldo 

Corporate  Research 
Operation  Department  (I) 

Chance  Vought  Aircraft,  Inc. 
Research  Division 
Box  5907 

Dal  las,  Texas  ( I ) 

Cross  Ma laker  Laboratories,  Inc. 

191  Mill  Lane 
Mountainside,  New  Jersey 
Attn:  Mr.  Richard  A.  Rossi  (1) 

Commanding  Officer 
U.  S.  Army  Signal  Res.  &  Dev.  Lab. 
Fort  Monmouth,  New  Jersey 
Attn:  SIGFM/EL-PRGd) 

Dartmouth  Co  1 1 ege 
Thayer  School  of  Engineering 
Hanover,  New  Hampshire 
Attn:  Prof.  R.  C.  Dean,  Jr.  (I) 

B.  F.  Goodrich  Company 
Research  Center 
Bucksvi I le,  Oh io 

Attn:  Mgr.  of  Res.  Operation  (I) 
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Goodr ich— High  Voltage  Astronautics 
Inc. 

Burlington,  Massachusetts 
Attn:  Dr.  C.  H.  Stockman 
Vice  President 

VIDYA  Associates 
2626  Hanover  Street 
Palo  Alto,  California 
Attn:  Mr.  D.  T.  Flood  (I) 

Warner  and  Swasey  Company 
Control  Instrument  Division 
32-16  Downing  Street 
Flushing,  Queens,  New  York 
Attn:  Mr.  G.  J.  Penzias  (!) 

Dr.  R.  Tourin  ( I) 

Prengle,  Dukler  and  Crump 
5417  Crawford  Street 
Houston  4,  Texas 
Attn:  Mr.  H.  Wm.  . /-angle,  Jr. 
Vice  President  (I) 

Aeronautical  Research  Laboratory 
USAF  Office  of  Aerospace  Research 
Wr ight-Patterson  A.F.  Base,  Ohio 
Attn:  Mr.  Erich  Soehngen  (2) 

Allied  Research  Associates,  Inc. 

43  Leon  Street 
Boston,  Massachusetts 
Attn:  Alan  Mironer 

Physics  Department  (I) 

Dr.  M.  P.  Freeman 
American  Cyanamid  Company 
1937  W.  Main  Street 
Stamford,  Conn.  (I) 

Engineer  ng  Projects  Laboratory 

School  of  Engineering 

MIT 

Cambridge  39,  Mass. 

Attn:  Gail  E.  Portridge  (I) 

Mr.  F.  C.  Gray 

Space  Technology  Laboratories,  Inc. 
I  Space  Park 

Redondo  Beach,  California  (I) 


